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Radiative Properties of Pair-instability Supernova Explosions 
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ABSTRACT 

We present non-LTE time-dependent radiative-transfer simulations of pair-instability super- 
novae (PISNe) stemming from red-supergiant (RSG), blue-supergiant (BSG) and Wolf-Rayet 
(WR) star rotation-free progenitors born in the mass range 160-230 M0, at 10^^ Zq. Al- 
though subject to uncertainties in convection and stellar mass-loss rates, our initial condi- 
tions come from physically-consistent models that treat evolution from the main-sequence, 
the onset of the pair-production instability, and the explosion phase. With our set of in- 
put models characterized by large ^^Ni and ejecta masses, and large kinetic energies, we 
recover qualitatively the Type Il-Plateau, Il-peculiar, and Ib/c light-curve morphologies, al- 
though they have larger peak bolometric luminosities 10^ to 10^° L0) and a longer dura- 
tion (~ 200 d). We discuss the spectral properties for each model during the photospheric and 
nebular phases, including Balmer lines in II-P and Il-pec at early times, the dominance of lines 
from intermediate-mass-elements (IMEs) near the bolometric maximum, and the strengthen- 
ing of metal line blanketing thereafter. Having similar He-core properties, all models exhibit 
similar post-peak spectra that are strongly blanketed by Fe II and Fe 1 lines, characterized by 
red colors, and that arise from photospheres/ejecta with a temperature of < 4000 K. Com- 
bined with the modest line widths after bolometric peak, these properties contrast with those 
of known superluminous SNe suggesting that PISNe are yet to be discovered. Being reddish, 
PISNe will be difficult to observe at high redshift except when they stem from RSG explo- 
sions, in which case they could be used as metallicity probes and distance indicators. 

Key words: radiation hydrodynamics - stars: atmospheres - stars: supernovae - stars: evolu- 
tion - stars: supernovae: individual: 2007bi, 2006gy 



1 INTRODUCTION 

About a decade ago, numerical simulations of star formation re- 
vealed t hat the first generat ion of stars may have been super- 
massive l lBromm et a"l]|2002h . holding great promise to discover 
them, not as stars, but r ather through their catast rophic explosions. 
iBarkat et al] \\96% and lRakavv & Shavivl j 19671) proposed that for 
stars sufficiently massive the production of e~e"'" pairs from 7- 
ray annihilation would lead successively to an implosion, a ther- 
monuclear runaway, and an explosion leaving no remnant behind. 
These events are termed pair-instability supernovae (PISNe). In the 
coming decade, the opportunity to identify PISNe and study them 
in detail may be realized with new observat ories such a s LSST 
and JWST (iScannapieco et al. 2005; Humme l et al.ll2012l) . How- 
ever, there is considerable uncertainty and speculation concerning 
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the first generation of massive stellar objects, impacting our under- 
standing of their formation, their evolution, and their explosion. 

Recent simulations for the formation of the first stars has 
significantly altered the original picture. The initial mass func- 
tion seems no longer to be heavily biased towards super-massive 
stars. In stead, it stretches down to lower masses, even very low 
masses JCaffau et al. 2012), perhaps exhibiting a bimodal distri- 
bution IsTac y et al. 201 1, 2012; Clark et al. 201 1). Binaries should 
also be created Istacv et aliboiol) . although the binary fraction is 
largely unknown. Simulations including feedback, in particular in 
the form of UV radiation, suggest that the circumstellar disk feed- 
ing mass into the primordial star is evaporated soon after the onset 
of steady thermonuclear burning. This feedback leads to the trun- 
cation of mass accretion and sets an upper mass limit for the newly- 
formed star of ~ 30 M0 jHosokawa et al.l201 ih . Theoretically, pri- 
mordial star formation does not seem to differ significantly from 
present-day star formation. Paradoxically, it fails to explain the ex- 
istence in the Local Universe of super-massive stars like 77 Car 
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jPavidson & Humphrevslll997l: Iffillier et al.ll200lh, massive stars 



in the Galactic centre jpiger e t al."2002VMartins et al."2008'). or 
the spectacular R136 nursery {Crowther et al. 2010). Since super- 
massive stars in the Local Universe systematically belong to dense 
star clusters, they may form not in isolation but instead dynamically 
and chaotically from merger events (e. g., iBonnell & Bate! [20051 : 
IVanbeveren et alj2009l ; |Pan et alj|2012bl) . 

More generally, the nature of superluminous SNe remains per- 
plexing. The interaction between two separate shells, which con- 
verts kinetic energy to thermal and radiative energy, is one way 
to boost the luminosity of a SN, because standard SNe typically 
have a hundred times more kinetic energy than time-integrated lu- 
minosities. Another way of producing a brigh t display is by de- 
positing magnetar radiat i on into the SN ejecta jMaeda et al.ll2007l : 
ikasen & Bildstenll2010l : IWooslevll2010h . The energy lost by the 
magnetar leads to the spin down of the magnetar which eventu- 
ally quenches the magnetar radiation. For a dipole field, the spin- 
down time scale is ts-p ^ A.^B^^Pl^d, where Bis is the mag- 
netic strength in units of 10^^ G and Pio is the rotation period 
P in units of 10 ms. Depending on the magnetar parameters, its 
spin-down timescale may be comparable to the decay timescale 
of ^^Ni, making it an attractive substitute for models that require 
a larger-than-average production of ^*Ni for powering luminous 
SNe. While a PISN may syn thesize a maximum of 57 M0 of 
^''Ni jHeger & Woosievll2002l) , it is unclear how much a "stan- 
dard" core-collapse SN can produce. This reflects the poorly known 
core properties at the time of death, and the viability of the neu- 
trino and magneto-rotational explosion mechanisms. The collap- 
sar model is one instance where copious amounts of ''''Ni may be 
produced ( IWoosle\][T993h , but the required core compactness to 
form such collapsars seem far more rarely matched by stellar evo- 
lution models, favori ng instead the formation of proto-magnetars 
jPessart et al]|2012dl) . 

Restricting our focus to very low metallicity environments the 
ansatz of the current study is that super-massive stars do form, pop- 
ulating the mass range 140-260 M©. Provided they do not lose 
much mass during the hydrogen-burning stage, they form a he- 
lium core with a mass in the range 64-133 Mq that should experi- 
ence the pair-production instability and le ad to a thermonuclear ex- 
plosion that com pletely disrupts the star jHeger & Wooslevll2002l : 
IWaldmaj2003) . How often this occurs is completely speculative as 
it depends on the unknown formation rate of the progenitors as well 
as their evolution. A major component of this uncertainty is mass 
loss since it determines whether the star will explode as an extended 
blue-supergiant (BSG), as a very extended red-supergiant (RSG) 
star, or instead as a very compact Wolf-Rayet (WR) star. This, in 
turn, impacts the supernova (SN) light curve morphology and thus 
the potential detectability and reliable identifi cation of these events 
jScannapieco et al1l2005l : iKasen et al.ll201lh . As we discuss here, 
a H-rich PISN from a RSG star at 1O~'*Z0 leaves little ambigu- 
ity about the origin of the event while connecting a superluminous 
Type Ic event with a H-deficient PISN is more challenging. It is 
unclear today if we have seen either H-rich or H-deficient PISNe 
jPessart et al]l2012ch . 

Mass-lo s s in m etal poor massive stars is highly uncertain. 
iBaraffe et al.l ( l200lh found that metal-poor massive stars are pul- 
sationally stable and, based on this analysis, should die with their 
hydrogen envelope. The metallicity dependence of radiation driv- 
ing in massive-star winds suggests that this other form of mass 
loss, critical for OB and WR stars at so lar metallicity, sho uld be 
largely inhibited at very low metallicity jCastor et al.lll975l) . As is 
quite typical for proto-stellar collapse in star formation, the col- 



lapsed c ore is endowed w ith a large amount of angular momentum 
at birth IStacv et al.'201 1"), which may affect its evolut ion in a num- 
ber of ways (Hirschi et al..,.2004a: .Hirschi, 2007: Heger & Wooslevl 



l201(]| : IChatzopoulos & Wheelej2012l : lYoon et al.l2012h . For exam 
pie, by increasing the oxygen core mass, fast rotation can lower 
the minimum main-se quence mass needed to encoun ter the pair- 
production instability jChatzoDoulos & Wheeleil2012h , thereby af- 
fecting the ^"^Ni-ma ss to ejecta-mass rati o and the resulting PISN 
radiation properties jPessart et al.ll20123) . High lumino sity poses 



a sev ere strain on the hydrostatic equilibrium of stars (Jo ss et al 



1971 , especially near break-up rotation speeds (JMaeder & Mevnetl 



2000l) . establishing ideal conditions for cent rifugally-driven mas s 



loss (a disk) and/or a continuum-driven wind jOwocki et al.ll2004l) . 
In this context, radiation- and/or rotation-driven mass loss may in 
fact compete with the episodic pulsations expect ed to occur in stars 
with a helium-core mass in t he range 40-63 M0 dHeger & Wooslevl 
l2002l : IWooslev et al.ll2007al) . 

Currently, the most promising PIS N candidate is SN2007bi 
jGal-Yam et al.ll2009l : lYoung et al.ll2010h . although this association 
raises several issues. The observations capture only the light-curve 
peak and beyond, making the explosion time uncertain, compro- 
mising the inference of the ejecta mass and ^^^Ni mass. Further- 
more, SN2007bi took place in an environment of 0.2-0.3 Zq; a 
hydrogen-deficient PISN at such a high metallicity seems d ifficult 
to accommodate with stell ar evolution ([Langer et al.| |2007|)FI The 
nebular spectral analysis of iGal-Yam et alj ( 20091) suggests a large 
ejecta mass, i.e., much larger than in standard core-collapse SNe, 
but the accuracy of this inferrence does not unambiguously support 
a PISN — within the errors these masses are compatible with a pro- 
genitor that would not experience the pair-pro duction instability. 

The alternative core-collapse scenario o f lMoriva et all ( l2010t) 
for SN2007bi requires extreme properties, which may be accept- 
able given the scarcity of these events, but it is also in contradiction 
with the environmental metallicity. Their proposed model yields 
a 36B (IB = IBethe — 10^^ erg) ejecta kinetic energy, which, 
combined with the binding energy of the progenitor, requires the 
remarkable extraction of < 50% of the gravitational energy of a 

typical neutron star. 

Kasen & B ildsteiJ ( 1201 Oh and lKasen etal] ( 1201 ll) have fueled 
this controversy further by obtaining equally satisfactory fits to the 
SN 2007bi bolometric light curve using either a magnetar model or 
a PISN model. In contrast, i n a recent study based primarily on the 
simulations presented here ( jPessart et al.ll20123) . we emphasized 
that PISN models systematically have red colors and narrow lines 
after light-curve peak - properties t hat are in strong cont radiction 
with the observations of SN2007bi. |Pessart et al. I ( l2012cl) suggest 
that delayed energy injection, as from a magnetar, into a lower mass 
ejecta would be more amenable to producing both blue colors and 
broad emission lines typical of such superluminous SNe. 

An independent problem with PISNe, if they were to be the 
first chemical nurserie s in the Universe, is that their distinctive 
nucleosynthetic yields jHeger & Wooslevl2002l : IUmeda & Nomotol 



l2002l : IChieffi & Limongi l2004ft are in apparent contradiction with 
those inferred from extremely metal poor stars dCavrel et al.ll2004l : 
ICohen et al.ll2008h . Such abundance studies suggest that the metals 



^ The lack of spatial resolution prevents the accurate inference of the 
metallicity at the actual site of SN2007bi, so one cannot exclude the 
possibility that the actual primordial cloud from which the progenitor of 
SN 2007bi formed may be of lower metallicity than the infeiTed value of 
0.2-0.3 Zq. In that case, a Type Ic PISN may not be excluded. 
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arose fro m "standard" core-collap se explosion of stars of moder- 
ate mass dHeger & WooslevlbOlOl) . although a contrived scenario 
involving a combination of S N with strong fallbac k and highly- 
energetic SNe was proposed bv lNomoto etal] ( |2006|) . 

In this work, we perform a quantitative study of PISN explo- 
sions arising from RSG, BSG, and WR star progenitors, investi- 
gating their photometric and spectroscopic sig natures, and docu- 
me nting in more detai ls the results presented bv lKasen et alj|20llh 
and lPan et alj ( l2012ah . We allow for non-LTE at all times and in- 
clude non-thermal processes, w hich are not just important for ex- 
citing helium atoms in SNe lb i lLucvlll99lh - they directly affect 
all species by driving level populations away from LTE. They also 
indirectly affect the thermodynamic state of the gas by channel- 
ing a fraction of the decay energy into excitation and ionization of 
ions and atoms rather than into heat. We compute the full time- 
dependent radiative transfer within the non-LTE framework, allow- 
ing the computation of the spectral evolution from the photospheric 
to the nebular phase. Multi-band light curves are obtained by inte- 
grating the resulting spectra over specified bandpasses. This unique 
approach allows to compute light curves and spectra simultane- 
ously and with the same level of complexity. It also allows us to 
start from physically-consistent ejecta resulting from stellar evolu- 
tion and explosion modelling. This detailed study aims at providing 
the key spectroscopic signatures of the three possible types of PISN 
explosions, with the hope of identifying the critical signatures that 
would help lifting the ambiguities surrounding the asso ciation of 
PISNe with events like SN 2007bi jPessart et alj2012cb . 

In the next section, we review the numerical setup for our cal- 
culations, including the pre-SN evolution from the main sequence 
(Sect. 12. It . the pair-instability explosion mechanism (Sect. 12. 2b . 
and the setup for the radiative-transfer calculations presented in this 
work (Sect. 12.3V We then describe the bolometric light curve of our 
PISN models, and compare them with those arising from "stan- 
dard" core-collapse SN explosions of 15-25 Mg, RSG, BSG, and 
WR stars evolved as single or binary stars (Sect.O. In Sect. |4] we 
coimect these light-curve properties to those at the photosphere. We 
then describe in detail the spectral evolution for each group of pro- 
genitors during the PISN photospheric and nebular phases (Sect.O. 
We compare our simulations to PISN candidates SN2007bi and 
SN2006gy in Sect. |6j additional information, in particular the 
proposition that the SN 2007bi li ght curve is more com patible with 
magnetar power, can be found in lDessart et al.l j20I2cl) . We finally 
advocate the potential use of PISN explosions as a diagnostic of 
the environmental metallicity in Sect. [7] We summarize our results 
and discuss their implications in Sect. [8] To limit the length of the 
main body of the paper, we provide additional results in an ap- 
pendix (provided as online material only), covering in more detail 
the evolution of the ejecta properties (e.g., ionization state and op- 
tical depth; Appendix A), the atomic data and model atoms (Ap- 
pendix B), and the dependency of our results on model atoms (Ap- 
pendix C). We finally provide tabulated values for the luminosity 
and magnitudes for all time sequences (Appendix D). 



2 PRE-SN EVOLUTION, EXPLOSION, AND MODEL 
SETUP 

The work presented in this paper was produced in several in- 
dependent steps. First, a large grid of massive-star progenitors 
was evolved from the main sequence until all models under- 
went collapse due to the pressure deficit following pair produc- 
tion. These simulations were then remapped into a ID radiation- 



hydrodynamical code, with allowance for (explosive) nuclear burn- 
ing and radiation transport, and followed through the explosion 
phase for a few years. A detailed discussion of the results from 
these two steps will be presented in a separate paper (Waldman 
et al., in prep.). Here, we focus on a subset of four PISN models 
from this large grid of models, and specifically arising from the ex- 
plosion of a RSG star, two BSG stars, and a bare He core. This 
paper discusses the evolution of these ejecta and their radiation 
taking into account both non-LTE and time-dependent effects. In 
Sect. l2. 1142.21 we present briefly the first two steps and then present 
the setup for the radiative-transfer calculations in Sect. 12.31 Addi- 
tional information is provided in the Appendix. 

2.1 pre-SN evolution 

Using the 1-D stellar-evolution code MESA jPaxton et al. Ilmib, 
we perform calculations of super-massive stars from the main- 
sequence until the onset of the pair-production instability. We fo- 
cus on progenitors with main-sequence masses between 160 and 
230 M0, adopt a metallicity of W^'^Zq and neglect rotation. In 
this paper, we focus on a small sample of PISN progenitors that 
synthesize a large amount of ^''Ni suitable to produce a superlu- 
minous SN. We thus ignore the less massive progenitors that en- 
counter the pair-production instability but produce low expl osion 
energies with a small ^''Ni mass (see, e.g.. lKasen et al.|[20I ih . We 
treat convection by the well known mixing length theory (MLT) 
using the Schwarzschild criterion and a mixing length parameter 
of 1.6, and ignore convective overshoot. During core hydrogen and 
helium burning the "basic" 8-species nuclear reaction network of 
MESA is used, including the following isotopes: ^H, "^'^He, ^'^C, 
"N, ^«0, ^'^Ne, ^^Mg. From the onset of core carbon burning, the 
nuclear reaction network is expanded. For the HelOO model (see 
details below) the a-chain elements are added, while for the other 
models, in order to account for more nuclei that may manifest in the 
resulting spectra, we use a more elaborate nuclear reaction network 
provided with MESA. In addition to the a-chain elements, it in- 
cludes the intermediate elements linking those through (a, p) (p, 7) 
reactions, namely: n, p, ^'"'He, ^^C, "'"N, ^''O, ^^F, ^°Ne, ^^Na, 

2*Mg, ^^Al, 28Si, 31p, 32s^ 35^^^ 36 39j^^ 40(^^^ 43g^^ 44^^^ 

*^V, ^**Cr, '^iMn, 52.54.56pe^ 55,56^^,^ 56 j^-^ Pair creation is treated 
in the eq uation of state used in MESA, as well as in Id (see below 
and, e.g., ICox&Giuiilll968h . 

A crucial ingredient in the models is mass-loss. Apart from 
the HelOO model, which was evolved without mass loss, we use 
two standard mass-loss prescriptions that are functions of the ef- 
fective temperature, Teft, of the stellar model. For cool stars, de- 
fined by Tcff < 10* K, we use the formula of Ide Jager et al.l 
1988h. multiplied by the metallicity-dependent factor [Z/Zq) 



Kudritzki et alJ 



19871) . For hotter stars, we use the formula of 



Ivink et aT 1 2001 ). which contains a {Z/Zq)°' metallicity depen- 
dence, where 0.64 < a < 0.69 depending on the effective tem- 
perature. Using these mass-loss prescriptions, the pair-instability 
is encountered when the star is a RSG. Our representative model 
for this scenario is R190. None of our models, evolved from the 
main sequence, lose enough mass to become a WR star, defined by 
-'fsurfaco < 0-4 and Tcff > 10'' K. This is expected - at the low 
metallicity of 10~* Zq adopted here, radiation-driven winds will 
be significantly inhibited during the object's life as an O star. 

To allow for uncertainties in mass loss (as well as for giant- 
eruptions, such as that associated with rj Car, which are not cov- 
ered by the above procedures), we perform additional simulations. 
To mimic the possibility of a stronger mass loss, we artificially 
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Table 1. Summary of the composition, in solar masses, for each model im- 
mediately before the onset of the pair-instability explosion. 



Model 


Mi 


Mf 


iH 


^He 


12c 


16o 


20Ne 


24Mg 


27A1 


28 Si 


R190 


190.0 


164.1 


23.6 


46.1 


5.3 


78.4 


8.6 


2.2 


0.06 


0.05 


B190 


190.0 


133.9 


5.7 


33.8 


5.3 


78.4 


8.6 


2.2 


0.06 


0.05 


B210 


210.0 


146.7 


3.7 


30.8 


5.9 


92.5 


10.8 


3.0 


0.08 


0.06 


He 100 


100.0 


100.0 


0.0 


8.5 


5.6 


75.4 


8.1 


2.1 


0.00 


0.07 



truncate the hydrogen envelope of the RSG progenitor immedi- 
ately before explosion, reducing the total mass and surface radius, 
and producing a BSG progenitor (models B190 and B210). To 
mimic an even greater mass-loss rate by means not understood to- 
day (although mass-transfer in a binary may be an option), we also 
evolve objects directly from a bare helium core (model HelOO). We 
summarize the properties for these four models in Table [T] In our 
nomenclature, RX, BX, and HeX refer to explosions as RSG, BSG, 
and He (WR) stars, where X is the initial mass. 

The evolution is generally followed using MESA's hydrostatic 
mode. When an hydrodynamic instability is encountered, which oc- 
curs invariably at central carbon exhaustion, the code is automati- 
cally switched to the hydrodynamic mode, and the run is continued 
until 2"Ne is exhausted at the stellar center. At this time, before 
^''O burning turns on, the profile is remapped into the ID radiation 
hydrodynamical code VID ijLivne 1993). We find that the nucle- 
osynthetic outcome is insensitive to the remapping epoch between 
hydrodynamic instability and neon exhaustion, but does change for 
later remapping epochs, mainly due to differences in the reaction 
networks used in the two codes. 

Our models typically reside on the main sequence for about 
2.5 Myr, gradually cooling from an O-star of about 80,000 K and 
~ 10 Rq, finally settling by the end of core helium burning, which 
lasts about 0.2-0.3 Myr, as a K-type supergiant of about 4300 K and 
4000-5000 Rq . Throughout its evolution, the luminosity of the star 
does not experience dramatic changes - it is 3~4x lO'' Lq, on the 
ZAMS, gradually growing by about a factor of 1.5 towards the end 
of the main sequence, and remaining virtually constant through the 
rest of the evolution. The star bums hydrogen and helium in shells 
for a period on the order of 10* yr, while its core contracts, until 
carbon is ignited at the center at Tc ^ Vf K, ~ 2x 10'* gcm~^. 
Since the carbon mass fraction is relatively low (on the order of 
0.06) the energy release cannot prevent the continued contraction 
of the core. Central carbon is exhausted in less than a year. An hy- 
drodynamic instability, which typically starts at Tc 1.2x 10^ K 
and pc ~ 6x lO^gcm"^, then occurs and this is followed by cen- 
tral neon exhaustion that takes several hours to be achieved. At this 
stage Tc ~ 2x 10^ K and pc ~ 4x 10^ g cm , while infalling ve- 
locities typically reach several 100 km s^*. The main properties of 
the models at the onset of the hydrodynamic instability are given in 
Table [T] We also show a Kippenhahn diagram for the RSG model 
R190 that illustrates the evolution of the internal structure and sur- 
face properties prior to the pair-production instability (Fig.[Tll. 

With our small grid, we thus cover the 3 different stellar types 
under which massive stars may explode, i.e., as RSG, BSG, or WR 
stars, resulting in explosions as SNe II-P, Il-pec, and Ib/c. The only 
configuration not considered here is SN Iln, which could stem from 
interaction bet ween consecutively eje cted shells following pair in- 
stability pulses jWooslev et al]|2007bl) . 



Table 2. Summaiy of the main properties for the four hydrodynamical input 
models B190, R190, B210, and HelOO explored in this study. All CMFGEN 
models are run assuming either local energy deposition, or non-local energy 
deposition (model name is then appended by "NL"). 



Model 


Type 


Mi 


Ml 




Ekin 


-"'-'cjecta 








[Mq] 


[M0] 


[R©] 


[B] 


[Mq] 


[Mq] 


B190(NL) 


BSG 


190.0 


133.9 


186 


34.5 


133.9 


2.99 


R190(NL) 


RSG 


190.0 


164.1 


4044 


33.2 


164.1 


2.63 


B210(NL) 


BSG 


210.0 


146.7 


146 


65.9 


146.7 


21.3 


HelOO(NL)" 


WNE 


100.0 


100.0 


1 


37.6 


100.0 


5.02 


HelOOK'' 


WNE 


100.0 


100.0 




40.9 


100.0 


5.00 



Notes: For all stellar-evolution models, the environmental metallicity is 
10^'* Zq. Type stands for the progenitor stellar type. The He-star model 
is run with two different sets of model atoms: HelOO uses the same model 
atoms as B190, R190, and B210. HelOOionl is identical in all ejecta and 
model parameters as HelOO except that it is modeled with the additional 
neutral species Mgl, Sit, S I, and Cat (Appendix C) . Corresponding 
properties of the HelOO model of lKasen et al.l )201 ll) , which we name 
HelOOK to avoid confusion. 

2.2 Explosion 

To model the explosion phase, we use a more extended reaction 
network, which in addition allows for (a, n) reactions up to cal- 
cium, for (n, 7) reactions up to nickel, and for the radioactive de- 
cay of ''^Ni -^^^ Co -s-'^® Fe. In total, we include the isotopes: n, 
p, *He, 12c, i*N, i«0 , 20 -22Ne, 23Na, 23-26Mg, 27ai, 2r-30si^ 

31p 31-34g SSqI 35-38^j. 39j^ ^^~*2q^ ''^Sc 47y 

Following central neon exhaustion, oxygen is ignited, and 
consumed in less than a minute. At this stage the typical central 
conditions are Tc ~3.5xlO^Kand/9c ~ lO'' gem , while the in- 
fall velocities reach several 1000 km s~i . In the subsequent ~ 10 s, 
central silicon is exhausted, and the implosion turns into an explo- 
sion. The central conditions reached at bounce depend on the mass 
of the model and the binding energy of the star, but are typically 
Tc ~ 3.5-6 X lO^Kandpc ~ 2-6x 10*^ gcm~^ 

During this explosive-burning phase most of the energy is re- 
leased by oxygen burning, and some by neon and carbon burning. 
For example, in our 190 IVIq model, ~ 40 IVIq of fuel are burnt, split 
between 30 IVIq of oxygen, 7 IVIq of neon, and 4 IVIq of carbon, 
with a total energy release of 44 B. This fuel is converted primarily 
into silicon (20 Mq) and sulphur (13 Mq), the remainder being ar- 
gon (3 Mq), calcium (2 Mq) together with ~ 3 Mq of ^^Ni. In con- 
trast, our heavier model B2I0 bums about 70% more fuel, namely 
50 Mq of oxygen, 9 Mq of neon, and 5 Mq of carbon, releasing in 
the process 75 B. This moderately increases the production of IMEs 
(23 Mq of silicon, 16 Mq of sulphur, 4Mq of argon, and 3Mq of 
calcium), but dramatically increases the ^®Ni production to 21 Mq. 

The energy liberated through explosive burning causes a huge 
increase in the ejecta internal energy (in the form of radiation) and 
kinetic energy. As the envelope expands and cools, radiative pres- 
sure gradients accelerate the material to its asymptotic velocity on 
timescales of hours to days. The energy lost to unbind the star is 
considerable, on the order of 10 B in all PISN models. Asymptoti- 
cally, the ejecta kinetic energy is 34.5 B for model B 190, 33.2 B for 
model R190, 65.9 B for model B210, and 37.6 B for model HelOO 
(Table [2l(. The mass- weighted average velocity for each model is, 
in the same order, 4200, 4000, 6100, and 5500 km s"* — this value 
can be used to infer a representative ejecta kinetic energy (that dif- 
fers from the exact value by <^ 30%) of 23, 26, 54, and 30 B. 
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Figure 1. Kippenhahn plot of the 190 Mq RSG model (i.e., R190), showing burning and convective regions, core sizes, mass, luminosity, radius and effective 
temperature as a function of the remaining time before the onset of the pair-production instability. The "luminosities" Ljj and Ljjo conespond to the nuclear- 
energy generation rates associated with H and He burning — there is also a strong neutrino luminosity (not shown here). Convective mixing refers to convective 
regions per se, while non-convective mixing corresponds to regions where semi-convection, overshoot and/or thermohaline mixing occur. 



Table 3. Summary of the chemical composition in ejecta models B190, R190, B210, and HelOO, immediately after explosive nucleosynthesis stops and prior 
to decay of unstable isotopes (specifically ^^Ni). We limit the table entries to H and He, CNO elements, the main IMEs, and finally Fe, Co, and Ni. The total 
mass of ^^Ni synthesized in the explosion is given in Table[2] The top half of the table shows the cumulative yields, and the bottom half the mass fractions at 
the progenitor surface. The ^"Ni surface mass fraction is zero for all models. All models assume an original metallicity of 10~* Z©. Numbers in parenthesis 
con'espond to powers of ten. 



Model A/h Mho Mc Mn Mq Mno Mms Msi Ms Mai Afca A/pc Mqo Mm 

[Mq] [Mq] [Mq] [Mq] [Mq] [Mq] [Mq] [Mq] [Mq] [Mq] [Mq] [Mq] [Mq] [Mq] 



R190 2.36(1) 4.53(1) 1.06(0) 5.06(-3) 4.81(1) 2.03(0) 2.87(0) 1.96(1) 1.30(1) 2.69(0) 2.36(0) 1.81(-1) 5.09(-3) 2.64(0) 

B190 5.72(0) 3.30(1) 1.05(0) 5.56(-3) 4.81(1) 1.99(0) 2.81(0) 1.96(1) 1.31(1) 2.72(0) 2.39(0) 2.03(-l) 2.53(-3) 3.01(0) 

B210 3.72(0) 3.04(1) 6.73(-l) 7.35(-3) 3.96(1) 1.43(0) 2.49(0) 2.34(1) 1.58(1) 3.52(0) 3.19(0) 4.35(-l) 6.15(-4) 2.15(1) 

HelOO 5.22(-7) 7.91(0) 1.36(0) 5.05(-3) 4.36(1) 2.11(0) 2.52(0) 1.94(1) 1.27(1) 2.54(0) 2.41(0) 2.67(-l) 6.27(-3) 5.04(0) 



Model Xh Xb_c ^Nc -''^Mg ^Si ^Ar -'^Ca ^Fc ^Co ^Ni 

R190 5.93(-l) 4.07(-l) 1.23(-7) 1.10(-6) 3.89(-7) 1.74(-7) 6.47(-8) 6.99(-8) 3.65(-8) 1.02(-8) 6.44(-9) 1.40(-7) 3.43(-10) 7.32(-9) 

B190 5.93(-l) 4.07(-l) 1.23(-7) 1.10(-6) 3.89(-7) 1.74(-7) 6.47(-8) 6.99(-8) 3.65(-8) 1.02(-8) 6.44(-9) 1.40(-7) 3.43(-10) 7.32(-9) 

B210 5.78(-l) 4.22(-l) 1.14(-7) 1.14(-6) 3.48(-7) 1.74(-7) 6.47(-8) 6.99(-8) 3.65(-8) 1.02(-8) 6.44(-9) 1.40(-7) 3.43(-10) 7.32(-9) 

HelOO 1.00(-15) 1.00(0) 3.37(-7) 9.87(-8) 9.16(-7) 1.74(-7) 2.76(-7) 6.99(-8) 3.65(-8) 1.02(-8) 6.44(-9) 1.41(-7) 3.43(-10) 7.32(-9) 
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Lagrangian Mass [M0] Lagrangian Mass [Mq] 



1.0 19.2 45.9 90.1 149.2 154.7 159.8 1.9 21.6 79.5 109.3 119.6 129.6 132.8 




Lagrangian Mass [M©] Lagrangian Mass [Mg] 

0.9 6.0 29.5 72.0 127.6 136.2 145.1 146.1 0.8 5.9 31.2 63.8 86.4 95.3 99.4 




Figure 2. Illustration of the chemical stratification in velocity space (bottom axis) and in mass space (top axis) for the four models discussed in this study 
(B190: Top left; R190: Top right; B210: Bottom left; HelOO: Bottom right). The stratification shown is not that of the original models - rather we illustrate 
the models used as input for CMFGEN in which the composition gradients were softened (except for B210). For better visibility, we limit the range and thus do 
not show the outer (homogeneous) regions. We also only show a selection of important species, rather than all species included in the calculations. 



The yields and the chemical stratification are described in Ta- 
ble [3] and Fig. |2] Ahhough PISNe are thermonuclear explosions, 
the bulk of the burning takes place in the inner envelope of the 
progenitor, yielding a chemical stratification more reminiscent of 
core-collapse Type Ic SNe: The ^^Ni is produced at the base of the 
ejecta, remains confined to the lowest expanding material with ve- 
locities less than 4000 km s~^ (Fig.[2ll, and the decay energy has to 
diffuse from there through the overlying mass to influence the outer 
ejecta and the light curve. 

The density profiles for each model, at the start of the 
radiative-transfer simulations, are shown in Fig. [3] We stitch an 
outer density region to all models to ensure an optically-thin 
outer boundary. With this procedure, the outer ejecta velocity 
for each model is 42000km s"^ (B190), 24000km s"^ (R190), 
50000km s"^ (B210), and 52000km s~^ (HelOO). The structure 
for model HelOO is quite smooth, with a density exponent of --^ 8 
beyond 5000 km s^^. There is a fair amount of structure in the 
density profiles for the BSG models and even more so in the RSG 
model, in part associated with the formation and propagation of 



reverse shocks at shell interfaces and exacerbated by the 1-D treat- 
ment. 

Hydrodynamical instabilities in PISNe have been investigated 
by|joggerst & Whalen ( 201 1). They find that mixing in Pop III 150- 
250 M0 SNe explosions is negligible compared to that obtained 
in 15-40 M0 Pop II counterparts, although their investigation does 
not contain H-deficient compact progenitors and they treat nuclear 
burning through an initial phase performed in ID. More recently, 
,Chen et al. (2012) have repeated such simulations and found that 
the treatment of nuclear burning in multi-D can enhance the level 
of mixing significantly. Since this issue is not settled, we decided 
to soften the composition discontinuities by mixing over adjacent 
mass shells, producing a mild smearing of the (high-resolution) ID 
input structures produced by Id. This is done more as a conve- 
nience than aimed at describing the physical effect of mixing seen 
in multi-D simulations. The chemical stratification of the models 
used as input to CMFGEN is shown in Fig.|2] Based on the present 
work, we however anticipate that mixing will not significantly al- 
ter the observables we describe here. Non-local energy deposition 
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Figure 3. Ejecta density distribution for the PISN Models B190, R190, 
B210, and HelOO at the onset of our CMFGEN calculations at 15.2, 33.4, 
15.3, and 10.5 d after explosion. For comparison, we oveiplot as a dotted 
line the slope for a power law with exponent 8 and 15. 



at very late times partly mimics the effect of mixing since 7-rays 
are then able to travel from their deeply located emission site in 
the ejecta to the 0-rich, the He-rich, and possibly the H-rich outer 
ejecta (Sect. [J!4l l. 



2.3 Radiative-transfer calculations 

Apart from the different hydrodynamic inputs, the numerical pro- 
cedure adopted in this work on PISN radiative-transfer simula- 
tions is the same as that employed previously f or the simulations 
of SN Il-peculiar an d in particular SN 1987A jPessart & Hillid 
I2OIOI: iLietalJlloTl), SNe Il-Plateau jPessart & HillieilllOllal) , 
SNe Ilb/Ib/Ic l lDessart et al.l201 ll , l2012bl) . lHillier & Dessarll l l2012h 
have recently given a full description of the code CMFGEN for SN 
applications, detailing the important components for the radiative 
transfer in SN ejecta. Since the code can naturally handle a wide 
range of composition mixtures (all important elements up to Nickel 
can be treated), the extension to treat PISN ejecta requires no spe- 
cific adjustment compared to the modeling of standard CCSNe or 
SNe la. 

We perform time-dependent simulations for the full ejecta 
at all times, and start by re-mapping the ejecta chemical stratifi- 
cation and structure (radius, density, temperature) computed with 
Id once homologous expansion is reached. In practice, we start 
the HelOO model at 10.5 d, the B190 model at 15.2d, the B210 
model at 15.3 d, and the R190 model at 33.4 d after explosion. This 
time increases with progenitor radius because for larger stars, the 
corresponding ejecta takes longer to significantly expand beyond 
its original spatial extent^ Given the long duration of the high- 
brightness phase for these events and the low early SN luminos- 
ity for the BSG and WR star explosions, starting this late is not a 



strong limitation here. For a description of the shock breakout and 
early l ight-curve signatures, the reader should consult [Kasen et al] 
l l201lh or Waldman et al. (in prep.)0 

While the ID Lagrangian hydrodynamical model computed 
with id contains about 2000 mass shells, the radiative transfer 
calculations are performed with a much smaller number of depth 
points. The critical scale to resolve is not mass but optical depth and 
we find that we obtain co nverged results with typi cally 5-7 points 
per optical-depth decade jPessart & Hilliedl20m) . In these PISN 
simulations, we use a fixed number of 150 depth points at all times, 
providing a satisfactory resolution of ^ 100 km s~^ near the base 
of the ejecta, decreasing to ~ 1000km s~^ at the fast-expanding 
outer edge. 

The bulk of the ejecta is located at high optical depth at early 
times and its base remains optically-thick until 200-300 d after 
bolometric maximum. However, for consistency and simplicity, we 
treat the entire ejecta in non-LTE at all times. This is not optimal 
in terms of memory requirements, but it turns out that the optically 
thick layers are close to LTE and converge very fast. It also allows 
us to smoothly handle the phase when the ejecta starts to thin out 
in the continuum - when this occurs some spectral regions may 
still be thick (e.g., in the UV where many lines overlap, blanket ef- 
ficiently the radiation, and drive the photons towards a blackbody 
distribution locally) while other spectral regions are already trans- 
parent (e.g., in the IR where the distribution of line opacity is more 
sparse). In general, numerous lines will remain optically thick for 
weeks after the continuum has become optically thin. Because we 
model the full ejecta, we adopt a zero-flux inner boundary. 

We explicitly treat non-thermal processes associated with ra- 
diatioactive decay. This is particularly important given the large 
mass of ^*Ni (we only consider the decay of ^®Ni and ^®Co in this 
study) and the low ionization conditions we obtai n in PISN ejecta. 
Through a solution of the Spencer-Fano equation jSpencer & Fand 
[1954I), we compute the contributions to non-thermal excitation and 
ionization for all elements as well as the contribution to heat- 
ing, and incorporate this in the non-LTE rate equations and the 
energy equation. A detailed prese ntation and app li cation of this 
non-the rmal treatment is given by iLi et al. I l l2012t) : I Dessart et al.l 
l l2012bl) . For the computation of non-thermal rates, we proceed by 
injecting all the decay energy as high-energy electrons of 2 keV and 
solve for the electron-degradation function at 2000 linearly-spaced 
energy bins down to 1 eV. At each of the 150 ejecta depth points, we 
construct and solve an upper diagonal matrix of 2000x2000, which 
is trivial, but does require up to 30 minutes in these PISN simu- 
lations. We reduce the computational burden of our non-thermal 
solver by considering in the computation only the most important 
ions of each species. As the electron density is the main quantity 
modulating the magnitude of non-thermal processes in a given cal- 
culation, and since it varies slowly between consecutive iterations, 
we recompute the non-thermal electron-degradation function every 
ten iterations. With this procedure, the treatment of non-thermal 
processes increases the total computational time of a model typi- 
cally by a few percent. 

The 7-ray escape from the ejecta is minute for all times con- 
sidered here, but the huge scale and mass of the ejecta combined 



^ In the RSG explosion model R190, the helium-rich shell, near 
4000km s^^, and sandwiched between the inner O-rich shell and the outer 
extended H-rich envelope has still not settled into homology. Here, enforc- 
ing homology leads to an artificial increase in velocity (and hence in radius) 
for the con'esponding layer, producing a higher ejecta mass by 7%. The CM- 
FGEN model R190 is then in fact 176 Mq. The BSG and He models are not 
affected by this issue. 



We note that given all the uncertainties on the PISN progenitors and the 
extreme sensitivity of the breakout signature on the progenitor atmospheric 
scale height, the presence of an optically-thick wind, or simply the stellar 
radius, the shock breakout signal is far from a clean signature to constrain 
the progenitor and explosion properties. 
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Figure 4. Synthetic bolometric light curves extracted from our non-LTE time-dependent CMFGEN simulations. Shown are our results for the PISN simulations 
discussed in this work: B190 (black), R190 (blue), B210 (turquoise), and HelOOionlNL (red). To better reveal the extraordinary properties of these PISN 
m odels, we a lso sho w our r ecent results for a "standard" SN II-P (olive; evolved at solar metallicity from a single 15 Mq main-sequence star; model sl5el2 
of IPessart & Hillier 201 1 j) . a "standard" Il-pec (orange; evolved at the LMC metallicity from a single 18 main-sequence star; model "lml8a7Ad" of 
iDessart & Hillieiti201(il) , a "standard" lib SN (green; evolved at solar metallicity from a binary 18Mq main-sequence star; model 365A2 of Dessart et al., 
in prep), and a SN la produced by the delayed-detonation in a Chandrasekhar-mass white dwarf with 0.67 Mq of ^^Ni (magenta; Dessart et al., in prep). 
In these last two models, 7-ray transport is solved for, predicting the fast decline rate during the nebular phase through 7-ray escape. Here, we use the term 
"standard" to stress that such II-P, Il-pec, and lib models match closely the bolometric light curve for the representative SN of each type observed routinely in 
the local Universe. Although B190(NL) and B210(NL) correspond spectroscopically to a SN Il-pec, R190(NL) to a SN II-P and HelOO(NL)/HelOOionI(NL) 
to a SN Ic, their photometric properties are obviously extraordinary. For each sequence, small filled circles indicate the actual post-explosion times at which 
the computations are performed. For reference, we give the bolometric magnitude on the right-hand-side ordinate axis. 



with the increasing 7-ray mean free path leads to non-local energy 
deposition within the ejecta. We find that this effect leads to the ap- 
pearance or strengthening of H I or He I lines at nebular times, and 
more generally tends to cause the broadening of line profiles with 
time as the SN progresses into the nebular phase. Hence, we typi- 
cally adopt local energy depos ition up to 200 d after explosion and 
employ a 7-ray transport code dHillier & D essart 2012;) to compute 
the non-local energy-deposition profile subsequently. The energy 
deposition from 7-rays is calculated once at the start of the simu- 
lation since it depends only on density, radius, time, and compo- 
sition (i.e., quantities that are not affected by the radiative-transfer 
solution or the ionization state of the gas) of the ejecta. Hence, 
CMFGEN models are run assuming either local energy deposition, 
or non-local energy deposition (model name is then appended by 
"NL", e.g., models R190 and R190NL). 

The model atoms ado pted for all simula t ions in this work are 
comparable to that used in iDessart & Hilliej ( I20T 1^ and are pre- 
sented in the appendix, in Table Bl, together with the relevant ref- 
erences. Not all simulations were performed at the same time and 
not all have the same model atoms. Simulations B 190, R190, B210, 



and He 100 were performed first, ran with the same model atoms, 
and neglected the neutral ions Mgl, Si I, Si, and Cal. We added 
these in the second-generation simulation HelOOionl. We in fact 
use that simulation to compare to observations of the Type Ic SN 
2007bi, as well as to gauge the blanketing effects from these neutral 
species. As we discuss in Appendix C, they have a significant influ- 
ence on the spectra and colors, but hardly influence the bolometric 
light curve. Hence, most properties of model He 100 apply to model 
HelOOionl; when these differences are significant, we specifically 
address them. 

The completeness of the model atoms is always an issue in 
our approach since we cannot blindly incorporate millions of lines 
from an extended list. Our line opacity enters the computations 
through the specific treatment of atomic/ionic levels and so we have 
to compromise between completeness and tractability. Each of the 
four simulations we present here requires about 40 steps each tak- 
ing about 3 days on 12 cores. The mem o ry is mostly taken by th e 
matrix to invert jHillier & Milled 1 19981 ; iHillier & Des"sar3l2012h . 
which contains all the terms from the linearized statistical and ra- 
diative equilibrium equations. The total memory allocated is typi- 
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cally iVT^ x (A^bands + 1) x ND where NT is the total number 
of equations (variables), A^bands is 1 (3) for a diagonal (tridiago- 
nal) operator, and ND is the number of depth points. We typically 
use iVr = 2000 (as in Table Bl), iVbands = 1, and ND = 150, 
corresponding to a memory of 2000^ x 2 x 150 x 8 ~ 9.6 Gb. 

Radiat ive-transfer simulat i ons fo r PISN e have been performed 
in the past. IScannapieco et all ( l200l l2005h focused on RSG star 
progenitors giving rise to SNe II-P, with approximations for the 
transfer such as flux-limited diffusion and blackbody spectra. Their 
bolometric light curves s hould therefore be accurate but their spec- 
tral colors are uncertain. lOal-Yam et alj 1^009) performed similar 
simulations but focusing on a range of He cores and explosion char- 
acteristics to model the observations of SN 2007bi (they model 
the light curve and o ne nebular-phase spectrum). As we discuss in 
iDessart et alj l l20123) , non-LTE line-blanketed simulations of PISN 
ejecta (which we present in detail in this work) yield a suitable 
match to the SN 2007bi light curve, but are in confl ict with nu- 
merou s observed spectral characteristics of SN 2007bi. [Kasen et al] 

l|) did an extended study of PISN light curves, including RSG, 
BSG, and WR star progenitors. As in our grid of models, they do 
not have a well motivated mechanism for predicting PISN explo- 
sions from BSG or WR stars so the final masses and stellar types at 
death remain rather speculative and thus represent exploratory ma- 
terial. With their radiative transfer, they considered both the multi- 
band light curves and spectra, highlighting for example the striking 
dependency on metallicity or the signatures for shock breakout, but 
their assumption of LTE for the gas, their approximate treatment 
of line opacity, and their focus on early-time photospheric-phase 
spectra limit the applicability of their models for confrontation to 
observations of PISNe candidates like SN 2007bi. This SN, like 
most superluminous events today, are still discovered at or beyond 
the peak of the LC, expected to occur weeks to months after explo- 
sion. At this phase, the SN ejecta has considerably expanded and 
starts becoming optically thin in the continuum, while the poten- 
tial presence of a large amount of IGE may maintain a high optical 
depth in lines for many months still. This configuration requires a 
detailed non-LTE line-blanketed treatment, which we adopt here. 

The radiative-transfer calculations performed with CMFGEN 
yield the level populations, electron density, and temperature as a 
function of depth and time. Summing over level populations gives 
the ion and atom populations as well as the ionization state of the 
gas. Concurrently, the properties of the radiation field are computed 
as a function of depth, frequency, and angle. In the next sections, 
we describe both the gas and the radiation properties, in an attempt 
to identify the photospheric-phase and nebular-phase signatures of 
PISNe. In the next section, we start by discussing the bolometric 
light-curve properties. 



3 LIGHT CURVE EVOLUTION 

One important output of CMFGEN is the emergent spectrum in 
the observer's frame. By integrating this flux over frequency at 
each time step, we recover the bolometric light curve. We show 
the synthetic bolometric luminosity for PISN models B190, R190, 
B210, and HelOOionlNL in Fig. |4l The fundamental signatures of 
these PISN models is the very long duration of the high-brightness 
phase, typically of a few hundred days, reaching peak luminosi- 
ties on the order of 10^" Lq at 150-200 d after explosion. Specif- 
ically, we find a bolometric maximum at 196, 220, 177, and 156 d 
for models B190, R190, B210, and HelOOionlNL. In the same or- 
der, the bolometric maximum luminosities are 4.4 x 10®, 6.3x 10^, 



3.0x 10^°, and 8.8x 10^ Lq. Allowing for non-local energy depo- 
sition changes these values by a few percent over this timespan. In 
the lower mass model HelOOionlNL, 7-ray escape starts being vis- 
ible at ~ 500 d after explosion, introducing a slight downward tilt 
of the fading rate below 0.01 magd~^, while for all other models 
(with an ejecta mass in the range 130-160 Mq), full 7-ray trapping 
holds exactly. Note that PISN explos ions may not syste matically be 
super-energetic and superluminous jKasen et alj|201 ll ), but our set 
of models was selected to be so, with kinetic energies of a few tens 
of B and nucleosynthetic yields including a few M0 of ^^Ni. 

The morphology of these light curves is analogous to the 
SNe we observe in the local Universe; the Type Il-Plateau class 
(R190), arising from a RSG progenitor jpalk & Arn"et^ll977h , the 
Type Il-pec class ( BI90 and B210) arising from a BSG progenitor 
( lAmettetalJl989l) , and the Type Ib/c class (HelOO) arising from a 
WR st ar progenitor l lEnsman & Wooslev 1988; Dessart et al. 20 iT], 
l20I2bl) . 

The anomalously large initial radius and kinetic energy of 
model R190 yields a "plateau" brightness over 10 times more 
luminous than obtained in the SN II-P model resulting from a 
1.2 B 15 RSG pro genitor (olive curve; model taken from 
iDessart & Hilliej|2011ah . The 2.63 Mg of '^'^Ni originally made in 
the explosion causes a strong re-brightening at 200 d after explo- 
sion, following the recession of the photosphere to deeper layers 
and the outward migration of a heat wave powered by ^®Co decay. 
This hump is never seen in SNe II-P produced from lower mass 
RSG progenitors because they produce too little ^®Ni in their ex- 
plosion. Beyond 350 d, model RI90 radiates at the instantaneous 
rate of decay energy deposition, even when allowance for 7-ray es- 
cape is made. The ejecta remains optically thick to 7-rays through- 
out the evolution modeled here, so that 7-rays may escape from the 
original site of emission but are eventually trapped somewhere else 
in the ejectaQ 

The light curve morphology for model B190 and B210 is 
explained i n the same terms a s for standard SNe Il-pec like 
SN 1987A l lBlinnikov et alJl2000l ; lDe"ssart & Hillieill2010l ; lLi et all 
l20I2h . However, as for R190, the light curve is much brighter and 
broader than for a standard Il-pec model (orange; model evolved at 
the LMC metallicity from a rotating single ISM© main-sequence 
star; lDessart & HillieJl20IOl ; lLi et alj|2012h . The luminosity falloff 
after shock breakout is much more gradual and extended and the 
onset of re-brightening, due to heating by radioactive decay, oc- 
curs later at about 50 d after explosion. The bright peak is much 
bigger due to the larger mass, optical depth, and hence diffusion 
time of the corresponding ejecta. In model B2I0, which synthe- 
sized 21.3 Mq of ^^Ni, the peak luminosity is 10*** erg s"* — ~ 7 
times larger than in model B 190 which synthesized only ~ 3 Mq of 
^''Ni. As in model R190, the larger ejecta mass prevents the rapid 
luminosity evolution seen in lower mass Type II SNe. 

The light curve morphology for model HelOOionI (and 
HelOO; see Appendix C) is typical of a Type I CCSN, with a faint 
post-breakout plateau, followed by a steep rise due to decay heat- 
ing at depth. The strong re-brightening starts ~ 30 d after explosion 
in this model, thus significantly later than in S Ne Ilb/Ib/Ic models 
tensman & Wooslevll988l ; lDessart et alj201 ll) . Provided mixing is 

* The use of a coarser grid and the successive re-mappings in the R190 
time sequence introduce a 10% change in the total mass and the ^"^Ni mass 
as time progresses. This causes model R190 to have the same nebular lumi- 
nosity as model B190, although the initial amount for ''''Ni is respectively 
2.63 and 2.99 Mq. This resolution problem affects negligibly models B 190, 
B210, and HelOO. See footnote^ 
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Figure 5. Evolution of properties at the electron-scattering photosphere, and specifically the radius (black; in units of 10^'' cm), the velocity (blue; in units of 
10000 km s^^), the temperature (red; in units of 10000 K), and the overlying mass (normalized to the total ejecta mass; green) for model B190 (top left), R190 
(top right), B210 (bottom left), and HelOO (bottom right). The dots con'espond to the actual times computed with CMFGEN. All ejecta models are optically-thin 
to electron scatteiing at > 550 d (see Tables A1-A4). 



weak in these PISNe l ljoggerst & Whaleij|201 ll : but see lChen et"al] 
l2012h . the rise to peak sh ould not start much earlier than obtained 
here jPessart et al.l2012bh . As for models B 190 and B210, the peak 
of the light curve is broad and luminous and qualitatively analogous 
between all three models. Although B190 and B210 correspond to 
a SN Il-pec, R190 to a SN II-P, and HelOOionl/HelOO to a SN 
Ic, their photometric properties are obviously special and it makes 
sense to retain the PISN calling for these models. 



4 PHOTOSPHERIC PROPERTIES 

As can be seen in Fig. |5] the photosphericlf] temperature is rather 
low in all models at all times, although a higher temperature is 
found for the R190 model because it starts from a much more ex- 
tended configuration that dwarfs cooling from e xpansion (associ- 
ated w ith the PdV term in the energy equation; iDessart & Hillied 
l2011ah . Since the photosphere is rather cool, PISN ejecta achieve 
a tremendous luminosity by remaining optically thick out to large 
distances - the maximum photospheric radii in models R190 and 
B210 exceed 10^^ cm, while the maximum photospheric radii for 
B190 and HelOO reach 5 x 10^^ cm. These radii are typically ten 
times larger than predicted for SNe resulting fr om the explosion of 
lower-mass massive stars, whatever their type jPessart et al.ll20o3 ; 
iDessart & Hiliieill2010l 1201 lal: IPessart et al.ll201lh. 



Figure |4] comprises all the SN types studied so far with CMF- 
GEN, and thus includes PISNe, SNe II-P, SNe Il-pec, and SN Ib/c. 
To complete the set, we add results from on-going SNe la sim- 
ulations for a delayed detonation in a Chandrasekhar-mass white 
dwarf synthesizing 0.67 M0 of ^^Ni (Dessart et al., in prep). In 
terms of peak brightness, such a SN la rivals luminous PISNe but 
its much lower mass prevents a sustained luminosity beyond ~ 30 d 
after explosion. 



^ Here, we define the photosphere as the ejecta location where the inward 
integrated electron-scattering optical depth is 2/3. At early times, when 
there is little line opacity, most photons will originate (i.e., be created) be- 
low the "electron-scattering" photosphere. At later times, the ejecta is metal 
rich and will have considerable line opacity and so using the Rosseland- 
mean or the flux-mean opacity would yield a location for the photosphere 
that is further out. Hence, our discussion of the "photosphere" is inerely 
iUustrative. 
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Figure 6. Evolution of the composition for representative species at the electron-scattering photosphere for model B190 (top left), R190 (top right), B210 
(bottom left), and HelOO (bottom right). ^''Nio refers to the sum of the mass fractions for the ^''Ni-decay products (it thus corresponds to the initial ^''Ni 
mass fraction in the coiTesponding mass shell). Most simulations have a short photosphere residence time in the H-rich envelope (B190 and B210), or none at 
all for HelOO in which the compact and H-deficient progenitor produces a photosphere in the O-rich shell up to 200 d. In model R190, up to 250 d is spent in 
the H-rich extended progenitor envelope, but by the time the photosphere reaches its base, the He-iich and 0-iich shells are fully recombined and transparent 
so the photosphere moves directly into the Si/S/lGE-rich inner ejecta. 



Beyond 50 d after explosion, and thus well before the light- 
curve peak for all models, the photosphere is at the recombination 
temperature for the dominant ion (H, C, or O), which happen to 
have a similar ionization potential of ll-14eV, and is on the or- 
der of 5000-6000 K. We thus expect red colors and spectra domi- 
nated by neutral and once- and twice-ionized ions (Appendix A & 
Fig. Al). 

As the evolution proceeds the photosphere generally migrates 
to deeper layers - this is best seen by looking at the photospheric 
velocity, or the mass that lies above the photosphere (Fig. |5j- The 
speed of recession is connected to the ability of the SN to remain 
hot and ionized because free electrons can contribute significantly 
to the opacity at the photosphere through electron scattering. Model 
R190 shows a slow progression of the photosphere inwards ini- 
tially, which remains in the outer 1% of the ejecta (in mass coordi- 
nate), before showing a prompt recession at the end of the plateau 
phase at ~ 230 d when it enters the core (this phase could have 
been better resolved with our simulations). The photosphere then 
crosses 60% of the ejecta mass in about 10 d as the ejecta becomes 
thin, essentially crossing the entire O-rich shell. In models B190 
and HelOO, the recession to deeper layers is fast early on and slows 
down as the decay heating becomes effective, halting temporarily 



the photospheric recession. In B210 the heat wave is so fierce that it 
makes the photosphere move outward tiirough ~ 27 Mq of ejecta 
material between days 50 and 100 (Fig. [5j. As a consequence of 
this phenomenon the photospheric velocity does not vary mono- 
tonic ally. 

The migration of the photosphere to deeper layers leads to the 
complete probing of the ejecta, from the progenitor surface to its 
core, over a time span of about 300 d. Because of the large ejecta 
mass and the presence of shells of distinct composition, exacer- 
bated by the low-primordial metallicity, this migration reveals an 
unprecedented chemical evolution at the photosphere (Fig.|6j. 

The enormous energy of a few tens of B released through 
burning (Sect. |2^ is compensated in PISN explosions by the huge 
ejecta mass. This is largely independent of the mass loss prescrip- 
tion since the pair-production instability can only occu r if the core 
mass is greater than --^ 60 M© jHeger & Woosle\j|2002l) . No matter 
what, the ejecta mass will be huge, well above that of any CCSN 
progenitor identified or inferred so far. A representative ejecta ve- 
locity is Urop = y/2E/M ~ 3200a/E^7a^io km s"\ where Ei 
is the kinetic energy in units of B and Mio is the ejecta mass in 
units of 10 M0 . A typical value for a RSG or BSG explosion in the 
local Universe is iSi = 1 and Mio = 1, i.e., v^cp ^ 3200km s^^. 
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Figure 7. Multi-band light curves for our PISN simulations obtained by integrating at each epoch the synthetic spectra over the band pass of filters U, B, V, R, 
I, J, H, and K. The absolute magnitude is shown. Magnitudes plotted here are also logged into Tables D1-D4. In each sequence, we obtain slight kinks in the 
magnitudes (especially in bands U, J, H, and K, where little flux is emitted) that do not occur in bolometric luminosity. These are tracked down to changes in 
parameters at the corresponding time, e.g., as the medium turns optically thin (switch to nebular inner boundary), change in SL assignment (split the lower 5 
to splitting the lower 10 levels for all ions), or change in model atoms (e.g., switch to a big Fe I model atom). While these glitches are not physical, the ejecta 
relax to the new proper solution at the next time step. 



For our PISN models, the increase in explosion energy more than 
compensates for the increase in ejecta mass and we have « rep ~ 
5000 km s~\ comparable to SNe Ib/c jPessart et ai]|201lh from 
intermediate-mass binary stars, but a factor of two lower than in a 
Chandrasekhar-mass white-dwarf explosion. Although both PISN 
and SN la are thermonuclear explosions, the bulk of the burning 
takes place in the inner envelope in a PISN, yielding a chemical 
stratification more typical of core-collapse Type Ic SNe: The ''"Ni 
is produced at the base of the ejecta, remains confined to the lowest 
expanding material with velocities less than 4000 km s^^ (Fig.[2l(, 
and the decay energy has to diffuse from there through the overly- 
ing buffer of mass to influence the outer ejecta and the light curve. 

The various properties described above for the bolometric lu- 
minosity and the photospheric properties assist the interpretation 
of the multi-band light-curves presented in Fig.[7] All models reach 
absolute visual magnitudes on the order of —20 to —22 mag. Model 
RI90 exhibits a quasi-plateau in the V band, together with a hump 
at 250 d caused by decay heating. The higher photospheric temper- 
ature makes the spectral-energy distribution (SED) peak short-ward 
of the Balmer jump for about two months following explosion. This 
evolution is typical of RSG explosions, with the magnitudes drop- 
ping fast in the blue but leveling for a long time in the red. After 
the secondary bolometric peak, aU magnitudes become fainter with 



time, but continue to ebb faster in the blue. The evolution of colors 
for models B190, B210, and HelOO are very different from model 
R 190 at early times. The fading after shock emergence is more ob- 
vious, and so is the re-brightening as the heat wave reaches the 
receding photosphere. The colors become bluer on the rise to peak 
before fading again as the ejecta becomes thin. This fading is again 
faster in the blue. 

An important feature of these multi-band light curves is that 
the colors are systematically red after the peak of the light curve. 
Figure[8]illustrates this color evolution for model HelOOionI, which 
is representative of BSGAVR star progenitors models at all times — 
for RSG star progenitor models, this similarity holds for a few 100 d 
after explosion and beyond. Before the peak, all models from BSG 
and WR star explosions are getting bluer: This is because the pho- 
tosphere is getting hotter due to the heating from decay at depth but 
without the "reddening" effects associated with high metal content 
at the photosphere (the photosphere still lies outside of these metal 
rich layers before the peak). The switch from red to blue and back 
to red as the SN evolves from the re-brightening to the peak, and 
eventually fading into the nebular phase is a distinctive signature 
of ''^Ni powered PISNe. This trend is somewhat similar to SNe la, 
with an additional re-hardening of the spectrum as the ejecta turn 
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Figure 8. Same as Fig. [T] but now showing selected colours for model 
HelOOionlNL. An important feature of PISN model HelOO/HelOOionI, 
shared with PISN models R190, B190, or B210 presented here, is the red 
colour of their spectra, especially at and beyond peak. 



transparent. We will return to this issue when comparing our PISN 
models to the observations of SN 2007bi (Sect. [611 . 



5 SPECTRAL EVOLUTION: GENERAL PROPERTIES 

From the preceding sections, we see that the four PISN models 
B190, R190, B210, and HelOO show significant differences in their 
light curves (peak brightness, time to peak, presence or not of a 
post-breakout plateau), which are easily connected to the funda- 
mental differences between the progenitors. We also expect the 
early spectra to differ significantly between models because of the 
differences between the progenitor envelopes. However, despite 
the different mass loss histories these models have similar super- 
massive helium cores and systematically produce large amounts of 
IMEs and IGEs from O burning (Sect. [2T[l2.2t . Consequently, we 
expect the nebular phase properties to be rather similar due to the 
similarities of the progenitor cores. Thus we describe the spectral 
evolution during the photospheric phase for each progenitor type 
separately (RSG, BSG, WR) before grouping them to discuss their 
nebular-phase properties. 

All our PISN models have a very large oxygen mass (39.6 
to 48.1 M0). However, as a consequence of the large ionization 
energy (13.6 eV) and the cool photospheric temperatures, O lines 
generally only have a weak direct influence on the spectra. The only 
exceptions are [O I] 6303-6363 A, which is seen in nebular-phase 
spectra, and a few lines of O I which weakly influence the early 
photospheric spectra of HelOO. 

Initially all models present a true continuum in the optical - 
that is the observed spectra consist of (broad) lines superimposed 
on an underlying background which is formed via continuum (i.e., 
bound-free and free-free) processes. However, eventually the pho- 
tosphere recedes into layers where IMEs (primarily Si and S) dom- 
inate the composition of the ejecta. At this, and later times, the true 
continuum flux weakens and then disappears. However, before the 
true nebular phase, a pseudo continuum, not generated by true con- 
tinuum processes, continues to be seen. It is generated by the huge 
number of lines associated first with IMEs and later on with IGEs. 
In some cases this allows, for as long as 500-600 d, the formation of 
P-Cygni profiles associated with line transitions from ions present 
in overlying mass shells. 



5.1 Spectral evolution during the photospheric phase of 
PISN model R190 

The large progenitor radius (~4000Rq) and explosion energy 
(33.2 B) for model R190 yield a SN which is photometrically 
(Fig. |4j and spectroscopically (Fig. |9} similar to a superluminous 
SN II-P. At early times, the outer ejecta and photospheric layers 
remain fully ionized (Fig. [5j Table A2) producing a blue contin- 
uum (red curve in the montage of Fig. |9j. The relative high ve- 
locity and steep density profile of the ph otospheric layers (Fig. [3] l 
yield broad but we ak P-Cygni profiles jPessart & Hillietll2005bl : 
IPessart et alj|2003) . At such times, the photosphere reside in the 
outer 1% of the ejecta mass, which is essentially pure hydrogen and 
helium (the metallicity is lO^^Zo). Consequently, we find only 
H I B aimer lines and a weak He 1 5 875 A. T hese are in fact typical 
spectral signatures of early-time SNe II -P jPastorello et alj|2006l: 
lOuimbv et al]|2007l : iDessart et ai]|2008h . but there are some criti- 
cal differences. To facilitate the analysis and line identifications we 
overplot in Fig. |9] (right panel) the synthetic spectra computed by 
omitting a given species|j 

As the ejecta expands helium quickly becomes neutral - the 
inefficient mixing of ^^Ni from gr eater depths prevent any non- 
thermal effects at the photosphere dDessart et al.ll2012b| ) at these 
times in these PISN models. The adopted metallicity of 10^* Z0 
impacts all species other than hydrogen and helium so that CNO, 
IMEs, and IGEs have little (direct) influence on the spectrum - up 
to 200 d after explosion, we can only see (identify) H I Balmer lines. 
As the ejecta evolve, we do not see the appearance of Fe II lines, 
Ca II H&K or Ca II triplet at 8500 A that are usually associated with 
Type II-P spectra. 

The lack of metals in the Hydrogen-rich envelope completely 
quenches the effect of line blanketing. As can be seen in Fig.|9] the 
continuum (red curve) follows closely the total synthetic flux, ex- 
cept over spectral regions where the H I Balmer lines reside. How- 
ever, by 200 d, a pseudo-continuum is visible in the full spectrum 
- a large number of lines, with no noticeable signature in the full 
spectrum, are contributing to the total flux, even in regions appar- 
ently devoid of lines. At this time, the spectral formation region is 
receding into layers containing IMEs. However, the photosphere is 
so cold that the emergence of line blanketing on the very weak flux 
short ward of 4000 A is difficult to see. By 300 d after explosion, 
the spectrum has retained a similar color to that at 200 d. How- 
ever, the pseudo-continuum is no longer produced by genuine con- 
tinuum processes (the red curve is down at zero flux) but instead 
by line opacity from Fe II. By 400 d, the main ion contributing to 
this background line-opacity switches from Fe II to Fe I and pro- 
duces an even redder pseudo-continuum (red and green curves in 
Fig.[9l right panel). By 600 d after explosion, the spectrum has vis- 
ibly turned nebular. 

In Fig. [To] (left panel), we illustrate the contrast between the 
spectrum from the explosi on of a 15 Mp RSG prog enitor at solar 
metallicity (model sl5el2: lDessart & Hillieill2011al) and that from 
the R190 model at 1O~*Z0. The line profiles are broader in the 



^ The other option is to include only the selected species. However, this 
procedure becomes dicey when continuum processes weaken and/or when 
line opacity dominates — both effects occur at nebular times in PISN explo- 
sions. In this case, taking out a dominant species completely alters the spec- 
trum and the result is unusable. So, for consistency, we adopt the same pro- 
cedure at all epochs and show these illustrations with the selected species 
omitted. This inanipulation is done only for post-processing and analysis — 
all CMFGEN simulations are obviously performed with all species included. 
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Figure 9. Left: Montage of synthetic spectra for the RSG model R190NL, showing the full spectrum (black) as well as the contribution from continuum 
processes only (i.e., bound-free and free-free processes; red). For completeness, we include an inset for the R190NL bolometric light curve. The emergent 
flux, scaled for better visibility (see inset for the absolute bolometric flux scale), is computed with CMFGEN by simulating the time-dependent non-LTE 
transport for the full ejecta from 36 until 1000 d after explosions. The sequence thus covers from times when the ejecta is optical thick until times when it is 
optically thin and nebular-line emission is strong. Right: Same as left, but we now overplot the corresponding spectra obtained by neglecting an ion (see color 
coding and labels at top-right). We illustrate the effect for ions that have a large impact on synthetic spectra. 
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Figure 10. The left panel shows a comparison between the RSG PISN model R190 and the equivalent explosion from a lighter massive star at Zq (15 Mq 
sl5el2 SN II-P model), while the right panel shows a comparison between the model B190 and the I8M0 "ImlSaVAd" SN Il-pec model at the LMC 
metallicity. The spectral moiphology is similar, except that the PISN explosions at 10""* Zq show no metals lines (not even Na I D or the Ca II triplet generally 
seen in all Type II SNe) at times when their equivalent solar/LMC explosions show strong line blanketing throughout the UV and optical ranges. 
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PISN model, but the contrast with the low-mass SN II-P model 
is small due to the modest increase in E/M. Although the col- 
ors and H I Balmer lines appear similar, metals visibly affect the 
low-mass SN II-P spectrum in the UV (strong line blanketing) and 
optical (weaker blanketing, primarily due to Fe II, Ti II, and Sc II) 
but have little effect on the PISN spectrum. This makes H/3 easily 
recognized as a P-Cygni profile in the PISN spectrum while in the 
solar-metallicity spectrum, H/3 overlaps with numerous metal lines 
(e.g.. Fell), and appear s only as an absorp tion feature. A similar 
effect was discussed bv lKasen et al.l 1 I2OI ih . Such RSG explosions 
in the early Universe could be used to constrain the environmental 
metallicity (or set an upper limit on the metallicity in case no metal 
lines are seen), either by means of the reduced blanketing in the UV 
or the w eakness of metal lines (e .g. Fe, Ca) in the optical (see also 
Fig. 9 of lDessart & Hilliej2005bh . The absence of Ca II lines at the 
recombination epoch is a unique feature of model R190, never seen 
in SNe II-P models, and never yet observed in SNe II-P. 

The Doppler velocity at maximum absorption in the broad H I 
Balmer lines tracks closely the decreasing photospheric velocity, 
although weaker Balmer lines from higher up in the series may un- 
derestimate it (Fig. ll II left column). This apparently peculiar find- 
ing is well explained by the line f ormation process in the homol- 
ogously expanding ejecta of SNe jPessart & Hilliej|2005al) . Non- 
UTE and time-dependent effects also contribute to maintain a sig- 
nificant line optical depth at large velocit ies, produci ng persistent, 
strong and broad lines for Ha and H/3 dUtrobin & Chugaii,2()0^ ; 
iDessart & Hilliej2008h . 

5.2 Spectral evolution during the photospheric phase of 
PISN models B190 and B210 

The spectral evolution of model B 190 and B210 reflects the typical 
explosions from BSG stars, objects characterized here by a sur- 
face radius of 186 and 146 Rq (Figs. [72UT3] ). As for model RI90, 
the large explosion energies of 34.5 and 65.9 B, respectively, pro- 
duce high representative ejecta velocities of ~ 5 100 km s^^ and 
6800 km s"\ well in excess of the value of 2800 km s ^ inferred 
for SN I987A jArnett et al.ll989l) . At early times, the spectral evo- 
lution for these two models is qual i tatively s i milar to the SN Il-pec 
model of lDessart & Hilliej l l20IOl) : lLi et"^ ( I2OI2I) and the obser- 
vations of SN I987A. 

Up to the phase of re-brightening (~ 50 d), the optical spectra 
are comparable to those of model R 190 at the recombination epoch 
and are composed exclusively of H I Balmer lines. The lack of met- 
als in the primordial gas, and the absence of mixing into the outer 
ejecta layers by enriched material (arising from steady or explo- 
sive burning), prevents any line blanketing in the UV or the optical. 
However, the smaller initial radius causes a much stronger cooling 
due to expansion, which induces early recombination at the photo- 
sphere and a fast recession of the photosphere (Fig.[6l(. While it took 
~ 250 d to reach the helium core in model R190, the photosphere 
recedes to the base of the He-rich shell as early as ~ 50 d after ex- 
plosion in models B 190 and B2I0. As seen before for model RI90, 
this corresponds to the appearance of new lines (from O I or Ca II) 
and the onset of line-blanketing effects, at this time primarily from 
IMEs. 

The photosphere recession to the helium core quenches the 
emission in Ha and other Balmer lines - a consequence of the small 
0(0. 1) hydrogen mass fraction in the he-rich shell and the lack of 
non-local 7-ray energy deposition into the H-rich envelope at such 
early times. 

In model B210, the larger ^''Ni mass causes a reversal of the 



photosphere, which moves back into the helium shell, but in model 
BI90, decay heating merely stalls the recombination wave halfway 
through the ejecta (in mass), within the 0-rich shell. It is around 
the light curve peak and beyond that the decoupling layers for the 
radiation recede to the IGE-rich layers, causing the appearance of 
iron line blanketing at ~ 200 d after explosion. 

In the right panel of Fig.[TOl we show a comparison, at the re- 
combination epoch , for models B 190 at 30 d after explosion and the 
SN Il-pec model of lDessart & HillieJ 1 I2OIOI) at 15 d after explosion. 
As for the RI90/sl5eI2 comparison, the PISN models B190 and 
B210 have systematically broader Balmer lines and lack the strong 
lines from Call, Nal, Sell, and Fell, present in the "lml8a7Ad" 
model at LMC metallicity. At such epochs these theoretical BSG 
explosions are relatively faint, and would be more difficult to de- 
tect than the much brighter RI90 model. 

For models B 190/B2I0, the Doppler velocity at maximum ab- 
sorption in the H I Balmer line is initially large, and declines as the 
photosphere recedes inward through the H-rich envelope (Fig.ll II 
right column, shown only for model B2I0). The decline is much 
more rapid than in R190, since the time scale for the photosphere to 
pass through the lower-mass H-rich envelope is ~ 50 d, much less 
than the ~250d for model RI90. In models B190/B2I0, the Hi 
Balmer lines present a maximum absorption at a Doppler velocity 
that strongly overestimates (by a factor of ~ 2) the photospheric 
velocity at early times - the same effect is seen in SN 1987 A 
jPessart & H illier 2010l). However, as the model re-brightens and 
its photosphere leaves the H-rich envelope, the maximum absorp- 
tion of the H I Balmer lines track a unique Doppler velocity, which 
corresponds to the velocity at the base of the hydrogen envelope. 

5.3 Spectral evolution during the photospheric phase of 
PISN model HelOO 

In this section, we describe the spectral evolution of a PISN arising 
from a WR progenitor using model HelOOionI, which differs from 
model HelOO by the treatment of additional neutral ions found to 
be important (Mg I, Si I, S I, and Ca I). We discuss specific differ- 
ences between HelOO and HelOOionI in Appendix C, which are 
essentially confined to the colors in the first few weeks following 
explosion, and additional line features and blanketing. 

At the start of the simulation, the HelOOionI and HelOO mod- 
els are at 10.5 d after explosion and the photosphere has already 
receded to the base of the 8 He-rich outer shell (Fig. |6j at 
~ 9000 km s^^ (Fig.|5). The composition is still He-rich but the O, 
Ne, Mg and Ca are also abundant, with mass fraction for IMEs that 
are orders of magnitude larger than for a solar mixture at a metallic- 
ity of lO"** Zq. This explains the dominant role of line-blanketing 
from IMEs, including Si and S, from such early times up to the 
peak of the light curve. It also explains the early appearance of 
strong Ca II lines. 

This chemical stratification, although not as obvious as in the 
H-rich models discussed in the preceding sections, is reflected in 
the spectral evolution of model HelOOionI (Fig.lI4b. Despite rep- 
resenting ~ 8% of the total mass of the ejecta. He is only present 
in the outermost ejecta shells. By 10 d after explosion, the photo- 
sphere has already crossed these He-rich regions so that there is 
no visible signature of He during the remaining part of the pho- 
tospheric phase. This is co nsistent with the (0.1) mass fraction 
of He in the outer layers jDessart et al.ll20ll1) together with the 
negligible no n-thermal excitation/ ionization at such epochs in the 
outer ejecta toessartetal]|2012hl) . We find Si I to be an impor- 
tant source of blanketing at the cool photosphere (for up to about 
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Figure 11. Top Row: Doppler velocity of the maximum ab.sorption (filled dots) in Ha (black), H/3 (blue), and H7 (red) for model sequence R190NL (left) 
and B210NL (right). For comparison, we overplot the evolution of the photospheric velocity (black line). Bott om Row: Evolution of the Ha line profile in 
model R190NL (left) and B210NL (right). Notice the typical P-Cygni profile morphology, already discussed in iDessart & Hilliej i2005bh . with a maximum 
absorption that may underestimate or overestimate the photospheric velocity, and the blue shifted peak emission. At nebular times, the line profile evolves 
towards a flat topped emission with a P-Cygni profile absorption fixed at ~4000km s~^ (R190NL) and ~8000km s~^ (B210NL), which in each case 
corresponds to the velocity at the base of the hydrogen envelope (see Fig.|2j- 



50 d). The continuum flux testifies for the presence of important 
photo-ionization cross sections (e.g., from Mgl with an impor- 
tant edge at 3757 A) and the absence of strong line blanketing 
in the optical, since the continuum flux closely follows the total 
flux for up to a month (Fig. I14t. To be complete about these line 
identifications, we find Ol lines at 6157, 6454, 7002, 7254, 7777, 
7990, and 9260 A; Mg I fines at 5167-5172-5183 A; Mg II lines at 
4481 and 7890 A; Si I lines at 3905, 4102, 5500-5800 (numerous 



lines), 5948 A, 7250, 7289, 7410, 7742-7799, 8648, 9413, 10585, 
10694, 10827, 10870, 10885, 10982 A; Sin lines at 4129, 5056, 
5957-5978, 6239, 6347-6371, 7849 A; Si lines at 6052, 6750, 
8712.44,8874.48, 9228, 9421, 9650-9681,10457, 10633; and Can 
H&K, 7103, 7162, 8133, 8201, 8248, 8498, 8542, 8662 A. The 
strong Mgl 5167 A is a striking feature of He 100 model around 
the peak of the LC. 

The larger photospheric temperature on the rise to bolometric 
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Figure 12. Same as Fig.|9] but now for BSG r 

maximum (Fig. |5j is visible in the hardening of the spectra (see 
also Fig. |7]l and the temporary weakening of blanketing from neu- 
tral IMEs. During that phase, the photospheric velocity is around 
5000km s~^, and the optical line profiles appear relatively nar- 
row (i.e., compared to the RSG/BSG progenitor models discussed 
above), apart from more optically-thick lines like Call. Numer- 
ous lines have a velocity at maximum absorption that underesti- 
mates the photospheric v elocity (Fig. I15t . As noted before and in 
iDessart & Hillieij l l2005ij) . this is caused by a projection effect in 
homologous ejecta since for increasing impact parameters, the lo- 
cation of maximum absorption tends to follow an iso-density con- 
tour, intersecting planes of ever smaller projected velocity. During 
this brightening phase, the continuum synthetic flux is again a good 
diagnostic for this heating and increased blanketing from earlier 
times. 



By ~ 150 d after explosion (i.e., the light curve peak), the pho- 
tosphere has now receded halfway through the ejecta, into the lay- 
ers dominated by Si/S/IGEs, and its temperature goes through a 
maximum (of ~ 6000 K). From then onwards, the spectrum red- 
dens, the blanketing from IGEs strengthens, initially through the 
effects of Fe II and eventually Fe I. The spectrum starts again show- 
ing strong lines of low-ionization and neutral species (as it did at 
~ 50 d), but with the additional contribution from metal line blan- 
keting. We also see the appearance of [Call] at 7291-7323 A. As 
seen before for the B190 and B210 models, the spectra at and be- 
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B190NL with 190 Mq on the main sequence. 

yond light-curve peak are red, not blue, with the bulk of the flux 
emerging at 5000 A, and progressing to longer wavelength as time 
proceeds. Throughout this phase, no Hel line is visible, which 
again confirms the notion th at a Type Ic classifica tion is not a guar- 
antee of helium deficiency ( iDessart et al]|2012bh . The fundamen- 
tal feature of these PISN spectra are the ubiquitous and concomi- 
tant presence of lines from O, Mg, Si, and Ca, in a much more 
pronounced fashion than generally seen in Type Ic SNe. The low 
metallicity of the environment is not as obvious as for the R190 
model, but noticeable after inspection from the lack of Fe II blan- 
keting in the 5000 A region; it is clear that the blanketing occurs 
but it is not IGEs like Ti or Fe that cause it here but instead IMEs 
and in particular Si. 

To illustrate more vividly the blanketing caused by the nu- 
merous lines of metals, as well as its evolution with time, we 
show in Fig. [T6] the distribution of the flux versus wavelength 
and velocity (equivalently the depth in the ejecta). Overplotted, we 
draw the photospheric location obtained for various opacity sources 
(electron-scattering, Rosseland-mean opacity, and flux-mean opac- 
ity). These photospheric velocities/radii are systematically smaller 
than obtained when accounting for the additional opacity of lines 
(red curve). Moving from the light-curve peak to ~ 300 d, the blan- 
keting is even stronger in the UV and optical, and it continues to 
strengthen as the ejecta cools and recombines. This blanketing en- 
hances the likelihood of interaction with a line for a photon emit- 
ted short ward of <: 5000 A. Through fluorescence, this photon can 
give rise to multiple lower energy photons, which, being subject to 
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Figure 13. Same as Fig.|9] but now for BSG model B210NL with 210 Mq on the main sequence. Notice, in particular, the presence of He 1 10830 A at nebular 
times and the red colors at all post-peak epochs, despite the huge ''''^ Ni mass synthesized in this explosion. 



a smaller opacity, can escape. Alternatively, blanketing lengthens 
the photon-residence time at depth, and thus enhances the proba- 
bility for absorption. Because of the decrease in temperature out- 
wardsQ this would translate into emission but now of lower energy 
photons. Both effects tend to deplete the flux where blanketing is 
strong, "reddening" the spectrum. Indeed, where the blanketing is 
strong, the flux is weak. 

5.4 Nebular phase spectral evolution of all four PISN models 

The transition to the nebular phase is not sudden, but occurs grad- 
ually as the intensity of weaker lines become less and less strong 
compared to the lines normally associated with the nebular phase of 
SN spectra (e.g., [O I] 6303--6363 A, [Ca II] 7290-7290 A and the 
triplet lines of Ca II at 8500 A; and Ha when H is present). 

This outer ejecta, which contains few tens of solar masses of 
IMEs, O, He, and in some cases H, is like a nebula on top of a 
centrally illuminating source. This persistent incoming radiation is 
complemented by non-local deposition of decay-energy from 7- 
rays escaping the core emission sites, which occurs when the 7- 
ray mean-free-path becomes sizable compared to the ejecta extent. 
This can resurrect the outer ejecta which had earlier on become 

^ In the optically-thin regions, the temperature distribution can be quite 
complex in practice due to density and composition changes, as well as the 
effects of non-local energy deposition. 



completely transparent and invisible. Hence, in PISN models that 
were faint initially, the nebular phase represents a means of probing 
the outer ejecta again. 

The strongest nebular line is [Ca II] 7290-7290 A - it has a 
moderate critical density of ^ lO*" cm"'', is easily excited by colli- 
sions with electrons, and for most of the SN evolution discussed 
in this paper Ca^ is the dominant ionization stage of Calcium. 
[Ca II] 7290-7290 A appears as a single broad feature at late times 
in all models - in fact it appears before the complete transition 
to the nebular phase, at 250±50d. In models B190 and B210, it 
is even present for a few weeks near ~50d. In RI90NL (which 
is representative of the whole set), we find that emission is lim- 
ited to the regions below 4500km s"^ in that model (except for 
Hq), thus limited at all times to the progenitor He-core. It peaks at 
4000 km s"'^ at 200 d and systematically recedes with time down 
to ~ 1500km s"^ at 1000 d. In contrast to Type II-SNe at near 
solar metallicity, there is no Call emission/absorption from the 
H envelope. This occurs because of the Ca deficiency outside of 
the core where Ca has a mass fraction of only 6.44 x 10"^. The 
low metallicity of the PISN progenitors in our set quenches any 
outer-ejecta (i.e., the H or He envelope) emission/absorption from 
species other than H/He. We note at nebular times the overlap of 
the [CaIl]7300A doublet with Cal lines at 7148 and 7326A (see 
Appendix C). 

The Call triplet is initially absent during the photospheric 
phase, appears as a relatively narrow P-Cygni profile as the pho- 



© 2012 RAS, MNRAS Q0Q.[Tl[36l 




Figure 14. Same as Fig. |9] but now for 100 Mq He-star model HelOOionlNL (same PISN model as HelOO, but the radiative transfer is computed with 
allowance for Mg I, Si I, S I, and Ca I; we also allow for non-local energy deposition). 




Figure 16. Variation of the co-moving frame flux versus wavelength and velocity (or depth in the ejecta) for model HelOOionI at 152 (left) and 296 d after 
explosion (right). We also overplot the location of the photosphere that results when including all sources of opacity (jagged red line), electron scattering 
only (bottom red line), the Rosseland-mean opacity (blue line), or the flux-mean opacity (green line). Finally, for illustration, we show on the right side of 
each panel the local energy deposition from radioactive decay of ^®Ni nuclei (black; in erg s~^ cm~^; normalized to its maximum value which occurs in the 
innermost shell) and the frequency-integrated CMF flux versus velocity (green). While the ^®Ni distribution is a given of the hydrodynamical input and thus 
fixed, the distribution of the energy deposition extends farther outward as time progresses and as the photospheres recedes. 
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Figure 15. Same as Fig. 1111 but now showing the results for model 
HelOOionlNL. We focus on early times when these lines remain optically 
thick. 



the time span 200-500 d, Ha shows a strong P-Cygni profile, in 
particular with a strong and broad absorption component. We show 
a montage for the complete evolution of Ha in Fig. [TT] for mod- 
els R190NL (left column) and B210NL (right column). At nebular 
times, the Doppler velocity at maximum absorption remains con- 
stant (see also the top panel in Fig.ll It and equal to the velocity at 
the base of the H envelope. The emission component is flat topped 
and its half-width also matches closely the velocity at maximum 
absorption. The H-deficiency in the inner ejecta (Fig. ^ together 
with the non-local energy deposition from 7-rays explains these 
properties. Note that due to strong line overlap with metal lines, 
Hq is not easily seen in all RSG/BSG models. Only B210 with its 
enhanced heating shows it unambiguously. 

In Fig. [17] we show the time evolution of the formation region 
for some of these lines in models R190NL and HelOOionlNL. The 
emission regions are ultimately bounded by the abundance distri- 
bution of the species under consideration, which reflect the pre-SN 
evolution, the regions where explosive burning took place, and the 
additional smearing we apply to all models (except for B210) — 
see Fig. |2] for details. For example, [O I] 6303-6363 A is present 
in all models at very late times. The lower the model E/M, the 
lower the expansion velocity of the 0-rich shell. We obtain repre- 
sentative velocities of 4000km s"^ (R190 and B190), 6000km s"^ 
(B210), and 7000km s^^ (HelOO). Oxygen emission thus arises 
from regions where both the O mass fraction and the mass den- 
sity are large. This emission is also influenced by optical depth ef- 
fects. In model R190, O emission comes from regions with veloci- 
ties extending from 1400 to 4500km s~^, while in model HelOO O 
emission is confined to regions with velocities between 3500 and 
10,000km s^^. At early times, the emission is weighted towards 
the outer regions because the high opacity limits the emission from 
the inner regions. At later times emission tends to be weighted to- 
wards the inner regions where the power source (the decay of Co) 
resides. 

In Fig. [TS] we compare this segregation of line emission at 
nebular times in model R190. The regions of emission for the 
main lines evolve little with time and reflect directly the chemi- 
cal stratification, yielding broad H I lines, intermediate witdh O I, 
and narrow Call. This is in contrast to solar-metallicity CCSNe 
at nebular times in which a significant fraction of the nebular line 
flux stems from the envelope unaffected by any nuclear burning , 
and thus at the prim ordial composition jPessart & Hillieill201 lal : 
iMaguire et al.ll2012l) . Similarly, in nebular-phase spectra of Type 
II CCSNe, Ha g enerally shows a narrow profi le and requires mix- 
ing into the core ( IUtrobinll2007l J Li et al ]|2012h . while here the Ha 
profile stems entirely from the outer H-rich sh ell. Mixing in those 
PISN e, which is expected to be inefficient jjoggerst & Whalen 
l201lU ut see also lChen et ai]|2012l) , should not alter this result. 



tosphere recedes into the helium core, before it turns into a set of 
relatively narrow emission lines at very late times. In model B210, 
the larger expansion rate causes the individual emissions to overlap 
and form a unique emission feature — different ejecta kinetic en- 
ergies and stratification yield different velocities for this emission 
region but the principle holds the same in all. 

Interestingly, some lines appear or re-appear at the nebular 
phase following the revival from non-local energy deposition. The 
most striking example of this is the re-appearance and persistence 
of Ha in model B210NL beyond 200 d, after the line had been 
weakening since the start of the simulations at 15 d (Fig.ll3t. Over 



6 COMPARISON TO OBSERVATIONS OF PISN 
CANDIDATES 

6.1 SN 2007bi 

Given the numerous uncertainties affecting the modeling of PISNe, 
our goal is to identify critical signatures from our simulations that 
will allow PISN candidates (e.g., SN 2007bi) to be confirmed or 
dismissed as a PISN. An important result of the present study has 
already been presented in the broader context of superluminous 
SNe, where we argue that PISN expl osions have probably not yet 
been discovered ( jPessart et ahlboiicl) . 
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Figure 17. Illustration of the spatial distribution of emission for representative lines (for multiplets, we show the bluest component) of important ions in model 
R190NL (top row) and model HelOOionlNL (bottom row). The quantity ^{v) is defined such that the hne flux is proportional to J' § d(logt)) (Hilliei 198'3). 
The time coverage for each panel is adjusted to capture the time when the line is unambiguously present in the synthetic spectrum. We find that all lines are 



primarily segregated according to the chemical stratification. 
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Figure 18. Comparison of the formation regions for Ha, [Ol]6300A, [CaII]7291 A, and CaIl8542A at 500d (left) and lOOOd (right) after explosion for 
model R190NL. 



The observed light curve of SN 2007bi reveals a peak at 
Mr ~ —21.3 mag and a slo wly decreasing nebul ar flux com- 
patible with full 7-ray trapping dOal-Yam et alj|2009l) . Our model 
HelOO (or HelOOionI; allowance for non-local energy deposi- 
tion does not influence the results for this comparison) matches 
closely these light-curve properties, although our m odel is fainter 
by ^ 0. 8 mag in the R band (Fig. (19). As discussed in lDessart et al.l 
( l2012d) . the models He 110 or Hell5 (corresponding to similar 
models to HelOO but from 1 10 and 115 Helium cores) match 
better the luminosity of SN 2007bi at and beyond the peak. How- 
ever, their spectral properties are comparable to those of HelOO. 
Overall, the explosion and light-curve properties that we infer for 
model HelOO are in agreement with the 100 Mq Helium-core mod- 



els co mputed by other groups dOal-Yam et al. I l2009l : I Kasen et all 

Soul)- 

Matching the light curve is a necessary step to validate a 
model, but it is not a sufficient test. Light- curve de generacy is epito- 
miz ed by the good qual ity fit obtained bv lKasen"& Bild sten ( 2010|) 
and lKasen et al.l ( 1201 ih using either a 20 ejecta powered by a 
magnetar or a 100 M0 PISN powered by 5 of ^^Ni. Modelling 
multi-epoch spectra is the critical next step to validate a proposed 
model. 

For a RSG progenitor, the high-brightness phase of the PISN is 
during the plateau, when the spectra are crammed with H I lines. SN 
2007bi is a Type Ic SN so it is unlikely to result from the explosion 
of a RSG star, even one endowed with few solar masses of ^''Ni 
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like our model R190. One may conceive of a contrived situation in 
which the SN 2007bi detection caught the ^^Ni bump that follows 
the plateau, at which time the H I lines may be weak, but why would 
the SN not have been detected during the preceding 200 d when the 
SN was even brighter in the optical? 

For the BSG progenitor B190, the light curve shape is in 
agreement with the observations of SN2007bi. The presence of 
hydrogen in the outer ejecta cannot be inferred from spectra at 
and beyond peak in this model (Fig. I12t so it would be at least 
in principle compatible with the observations of SN 2007bi and its 
Type Ic classification. Model B210 does show Ha unambiguously 
at nebular times (Fig.|13tFlbut it produces a too high luminosity for 
SN2007bi. In model B190, the presence of an hydrogen envelope 
(not directly inferred from spectra at and beyond peak) leads to a 
smaller expansion rate of the He-core than in the HelOO model. The 
broad line features seen in the optical spectra of SN 2007bi suggest 
fast expansion of the emitting layers and so, in this PISN context, 
would favor the lowest mass progenitor capable of producing 3- 
5Mo of ^'^Ni. Hence, our model HelOO (HelOOionI) seems the 
most suitable of all four for this PISN candidate and we therefore 
focus on this below. 

The main spectral signatures for model HelOOionI are the 
dominance of absorption and emission processes from IMEs in the 
spectra up to the peak, and a hardening of the spectrum on the rise 
to peak (Fig. I14t . After the peak, the spectrum formation region 
is located within the Si/S/Fe-rich layers. The high metal content 
together with the cooling of the photosphere, leads to strong blan- 
keting in the blue and a severe reddening of the spectrum. This 
alone is the most funda mental disagreement with the observations 
of SN2007bi (Fig [HI lOessart et al] 1201 2d) , and in fact of other 
super luminous SNe which seem to be systematically blue after the 
peak dOuimbv etalj201ll) . Our PISN simulations are in contradic- 
tion with this observed property. Indeed, it would be surprising for 
PISN explosions to retain hot photospheres free of blanketing few 
hundred days after explosion. 

On the rise to peak, the photosphere feels the heat wave pow- 
ered by decay and diffusing from greater depth, and consequently 
cannot feel the blanketing effects of these buried IGF; the spec- 
tra before peak are bluer. After the peak, the photosphere cools and 
recedes to IGE-rich layers, feelin g the full effect of m etal-line blan- 
keting. The comparison made bv lYoung et al. I( l2010l) between typ- 
ical SN la spectra and those of SN2007bi reveals a striking sim- 
ilarity, but this similarity suggests SN 2007bi is not a PISN. In 
the best fitting model HelOO of a 100 Mq ejecta, the ^'^Ni mass 
is 5.02 M0, thus 5% of the total mass. This is closer to the value 
of --^ 1% characterizing SNe II-P t han the value of ~ 5 0% char- 
acterizing SNe la. As discussed in iDessart et al. I l l2012cl) . models 
with increasing ^^Ni-to-ejecta-mass ratios have bluer colors but 
even in the most extreme case of PISN model Hel25, this ratio 
is 0.26 and the colors are still too red to match those of SN2007bi 
— this Hel25 model is also much too bright with its ^''Ni mass 
of 32M0. In the models HelOO — Hell5, which match quite well 
the SN 2007bi light curve, the amount of energy released per unit 
mass is not so favorable to produce a hot emitting ejecta at and be- 
yond the peak of the light curve, as suggested by SN 2007bi. And 
indeed, our nebular-phase synthetic spectra are systematically cool 
and suffer severe line-blanketing dwarfing the flux short ward of 
~ 5000 A. This color mismatch can be inferred by comparing the 
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Figure 19. Top: R-band light curves for PISN models B190. B210. and 
HelOO, together with the observed R-band light curve of SN2007bi (ad- 
justed horizontally for convenience). The abscissa is the time since bolo- 
metric maximum — it differs by merely a few days from the time of R- 
band maximum used below. No correction fo r reddening is applied . Middle: 
Montage of observed spectra of SN2007bi iOal-Yam et"anl2009l) . stacked 
arbitrarily for visibility. Labels refer to the time since R-band maximum. 
Bottom: Montage of synthetic spectra for model HelOOionlNL at contem- 
poraneous epochs. 



° This results from the huge °Ni production combined with non-thermal 
ionization and excitation, as well as non-local 7-ray-energy deposition 



© 2012 RAS, MNRAS 00Q.[Tlt36l 



Simulations ofPISNe 23 



model He 100 color B — R with the observations of SN2007bi. In 
the model, it is ~ 0.55 at the light curve peak and subsequently, it 
steadily rises to reach 2.5-3 after 200 d (Fig.[8l(. In contrast, the ob- 
servations of SN 2 007bi show a- B-R of - 0.23 mag (0.9 mag) 54 d 
(150 d) after peak jYoung et alj201(t) . hence a much bluer colour. 

A further problem with the similarity between S Ne la spec- 
tra and SN 2007bi illustrated by lYoung et all ^20 id) is that in a 
PISN explosion, the ^''Ni is produced in the hottest and densest 
regions of the ejecta. In model HelOO, these regions are ejected 
at a speed 4000 km s~^, while in SNe la these layers reach up 
to 15000 km s~^. This strong contrast in expansion rate for the 
^^Co/^®Fe-rich layers conflicts with this spectral similarity: metal 
lines in SNe la are broad while those in PISNe are narrow. It is 
al so illustrated by the in consistency of the PISN model presented 
in lGal-Yam et alj l l2009l) : inferring a mass- weighted mean velocity 
of 8000km s~^ and an ejecta mass of 100 M©, the explosion has 
to deliver ~ 100 B (note that ~ lOB is needed to just unbind such 
a super-massive compact star), which is then incompatible with the 
inferred 3-10 M© of ^®Ni synthesized in the explosion. For com- 
parison, the extreme model B210 presented here releases "only" 
75 B through burning, producing an ejecta with a kinetic energy of 
65 B^ 

lKasenetalJj201lh have compared one synthetic spectrum for 
their HelOO model at 50 d before peak with the spectrum of SN 
2007bi at 54 d after peak. They do reproduce the blue spectrum 
but this stems from the fact that they have a 100-d mismatch. At 
50 d after explosion, the color is no longer blue and the spectra are 
strongly line blanketed, in contradiction with the observations. 

It thus seems unlikely that SN 2007bi is a PI SN explosion. In 
addition with the issues raise d by star formation jHosokawa et al.l 
I2OI ih and stellar evolution jLanger et al.l l2007h for such super- 
massive star progenitorslfl we have provided independent argu- 
ments that suggest that even if PISN ejecta existed, their radiative 
properties would conflict in nu merous ways with the observations 
of SN2007bi. As we propose in lOessart et al.l ( l2012d) . we instead 
favor a magnetar scenario for SN 2007bi, as well as for superlumi- 
nous SNe characterized by blue colors. 

6.2 SN2006gy 

The spectral evolution of our H-rich PISN models (i.e., R190, 
BI90, B210) can also be compared to that of PISN candidate 
SN2006gy, whose superluminous display was associated with 
^fNi-power or CSM interaction jOfek et all |2007| ; ISmith et al.l 
l2007h . This superluminous Type Iln SN reveals a spectrum dom- 
inated noticeably by Hi Bal mer and Fell line s with narrow and 
broad velocity components l lSmithetal.ll20IOl) . Its spectral mor- 
phology suggests a slow cooling of the emitting region, with little 
signs of line blanketing, and, surprisingly, no emission of the usual 
[Ca II] and [O I] nebular lines seen in massive-star explosions. 

The ^''Ni-power vs. CSM-interaction models can be discussed 
further in light of our simulations. The explosion of a RSG progen- 
itor is directly rejected from the bell-shape light-curve morphol- 
ogy of SN 2006gy, while the presence of hydrogen would require 
a blue or a yellow supergiant (in agree ment with the expe ctations 
from stellar evolutionary calculations; iLanger et al.ll2007l) . How- 
ever, our simulations, wherein the luminosity is powered by ra- 
dioactive decay primarily, indicate that by the time the SN reaches 

^ On possible way to overcome the evolution argument is throug h mergers 
in a dense star cluster, as recently discussed bv lPan et al.l i2012bl) . 



its peak brightness, the photosphere is in the Helium core and radi- 
ates a spectrum that contains little sign for the presence of hydro- 
gen (models BI90 and B210). The bulk of that energy arises from 
the core and the emission is channeled for a sizable part into those 
fine structure lines, which should thus be observable. The fact that 
SN2006gy retained a Type II spectral morphology for hundreds 
of days with no signs of blanketing short ward of 5000 A and no 
forbidden line emission from [Ca II] and [O I] invalidates the ^^Ni- 
power model for SN 2006gy. 



7 PISN AS METALLICITY INDICATORS IN THE 
UNIVERSE 

Being so luminous, PISN explosions represent attractive probes of 
the young Universe where the first stars formed. In particular, those 
that die as RSG stars would sustain a high luminosity in the op- 
tical from early after breakout until the transition to the nebular 
phase about a year later (in the SN frame). As discussed earlier 
(Sect. I5.lt , the representative model R190 has a nearly pure H I 
spectrum (with small contributions from He I) with no sign of line 
blanketing from either IMEs or IGEs for ^ 250 d after explosion. 
During that extended time, which may last a few years for an ob- 
server on earth depending on the redshift, spectroscopic observa- 
tions could reveal the environmental metallicity out of which the 
PISN progenitor formed. At large redshift, th is determination i s 
generally done from analysis of nebular lines l lOsterbro^ll989h . 
but doing this task using the SN spectrum offers a very interest- 
ing alternative. In the present case reported, the metallicity is so 
low that no metal lines are seen at early times, which would be an 
unambiguous signature that the metallicity is extremely low. For 
higher metal abundances, one ex pects a gradual rise in associated 
line strengths l lKasen et al.ll201ll) . 

As we discussed in detail, in particular for model HelOOionI, 
on the rise to peak the photosphere probes the He-rich and O-rich 
envelope above the IGE-rich inner ejecta and thus still reflects the 
primordial abundances for IGEs. After the peak, this information is 
all lost since we see emission/absorption from the ashes produced 
by steady and explosive burning. This still leaves an extended win- 
dow for determining iron abundances, for example, for the primor- 
dial gas out of which the star formed. 

In Fig. |20] we show how the iron mass fraction varies at the 
photosphere for all models, with horizontal lines indicating the lev- 
els corresponding to the SMC, the LMC, and the solar iron mass 
fraction. Over a few hundred days, the iron mass fraction varies in 
all cases over 7 orders of magnitude, i.e. from the iron-deficient 
outer ejecta to the pure iron layers of the inner ejecta. Note how- 
ever that the primordial metallicity is preserved at the photosphere 
of model R190 for ~200d (because the thermalization region is 
much deeper than the electron-scattering photosphere, metals in- 
fluence the emergent spectrum of model RI90 before 250 d; see 
Fig.H. 



8 CONCLUSION 

In this paper, we have explored the radiative properties of PISNe 
arising from the explosion of RSG, BSG, and WR star progenitors 
at a very low metallicity of 10"* Zq and without rotation. Our ap- 
proach is based on physically-consistent modeling of the evolution 
from the main sequence until the onset of the pair-production insta- 
bility using the stellar evolution code MESA, the explosion phase 
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Figure 20. Illustration of the iron mass fraction at ttie electron-scattering 
photosphere versus time since explosion for each model. Two competing 
effects cause this quantity to generally increase with time, i.e., the photo- 
sphere recession to IGE-rich layers as time proceeds and the decay with 
time of ^®Ni and ^^Co isotopes. We draw a dotted horizontal line to rep- 
resent the iron mass fraction in the SMC, LMC, and in the solar neighbor- 
hood. More generally, RSG-star explosions constitute an attractive probe to 
constrain the metallicity, especially those occurring in the early Universe. 



and the evolution of the ejecta with the radiation-hydrodynamics 
code id (with allowance for explosive burning), and the non-LTE 
time-dependent radiative-transfer modeling with CMFGEN for the 
multi-band light-curves and spectra from the photospheric to the 
nebular phase. In this paper, we have focused on this last step, while 
the former two steps will be described in Waldman et al. (in prep.). 
In practice, we focus on a 190 M© RSG star (model R190), a 190 
and 210 Mq BSG stars (models B190 and B210), and a 100 Mq 
WR star (model He 100). 

The exceptional progenitor and explosion properties of the 
PISNe we modeled set them apart from the SNe we routinely ob- 
serve in the local Universe. With ejecta masses of ^ 100 , ejecta 
kinetic energies of few tens of B, ^^Ni yields of a few Mq, their 
luminosity peaks at ^ 10^" about 150-200 d after explosion, re- 
maining above ~ 10^ for ^ 100 d. Their large ejecta mass en- 
sures full 7-ray-trapping so that the ^^Co decay rate gives the dis- 
tinctive ~ 0.0 1 mag d^ ^ nebular decline rate up to ~ 500 d after ex- 
plosion (~ 300 d after bolometric maximum). Over the next 500 d, 
the fading rate slowly steepens in the lower ejecta-mass model 
HelOO but stays the same in the more massive ejecta models. 

The ingredients ruling the evolution of "standard" SNe in the 
local Universe apply to PISNe directly. The influence of the progen- 
itor radius is key in yielding the high-luminosity plateau of the RSG 
model R190, while the more compact BSG and WR progenitors are 
characterized by a single highly luminous '''^Ni-powered peak. Fur- 
thermore, the large progenitor radius of model R190 yields a hot 
ionized photosphere for weeks after explosion. In the BSG model 
B190, recombination turns on much more quickly, as in SN 1987A. 
In the helium-rich model HelOO, the re-brightening phase to the 
bolometric peak corresponds to a photospheric residence in the O- 
rich core, with inefficient non-thermal excitation of the outer he- 
lium. Hence, our model R190 would be classified unambiguously 
as a SN II-P. Models B190/B210, despite having a si zable amount 
of H a nd He, would be Type lie if observed early on jPessart et al.l 
l20I2bl) but most probably a Ic if discovered at peak or beyond. 



Finally, model HelOO would definitely be a Type Ic despite the 
~ 8 M0 of helium in its outer ejecta. 

Despite these significant spectral and light-curve differences 
up to the bolometric peak, all models reveal very similar light- 
curve and spectral properties beyond it. This naturally stems from 
the comparable core properties of the progenitor stars — here, dif- 
ferences in cumulative mass for O or ^''Ni change little the tem- 
perature or ionization controlling the radiative properties from the 
core and conspire to produce this degeneracy. 

Tw o fundamental prop erties emerge from PISN simulations 
(see also lKasen et alj|201 ih . First, homologous expansion in such 
massive ejecta implies that while the outer layers may be travel- 
ing fast, as witnessed early on from the width of the Balmer lines 
in the SN II-P model R190, the inner ejecta travels at relatively 
small speeds and should therefore lead to rather modest nebular 
line widths. The ^''Ni rich core travels at < 2000 km s^^ in model 
B190/RI90 and at <: 4000 km s"^ in models B210/Hel00. Hence, 
although these objects are thermonuclear explosions, their chemi- 
cal stratification is closer to that of core-collapse SNe than SNe la. 
In practice, we find all our PISN simulations to follow this principle 
and reveal nebular lines with a narrow width, in particular arising 
from the 0-rich and Ca-rich shell on top of the '"''Ni-rich core. 

The second, perhaps counterintuitive, property of these PISNe 
is the relatively modest ejecta and photospheric temperatures of 
< 5000 K at and beyond the bolometric maximum. What controls 
the ejecta ionization and temperature is not just the mass of ^"^Ni, 
which is admittedly large, but instead the ratio of '''^Ni mass to 
ejecta mass, which is typically quite small in these PISNe, again 
more typical of core-collapse SNe than SNe la. Indeed, the large 
ejecta mass means that the energy released is to be shared. This 
is an inevitable property of PISNe since pair-production can only 
occur in stars with super-massive He/CO cores. Indirectly, this 
increases the diffusion time of the ejecta, which becomes trans- 
parent late, at < 200 d in this set. By then, the ejecta has ex- 
panded to lO'^'' cm, corresponding to an increase in radius of 10^ 
for the core material at that time. At such late times, ^®Co decay 
dominates, with a weaker heating rate compared to ^®Ni. Conse- 
quently, the PISN spectra at and beyond bolometric maximum are 
red, i.e., little flux is emitted shortward of ~ 5000 A, and strongly 
blanketed by Fell and eventually Fel. This is in stark contrast 
with any known sup erluminous SN today. As we emphasize in 
lOessart et alj l l2012cl) , these properties are a t odds with the obser - 
vations of the PISN candidate SN2007bi ( Gal-Yam et alj|2009h . 
An attractive alternative, which needs further study, is the magne- 
tar model whereby energy is added to a SN ejecta but without the 
compromising blanketing effects of metals. 

The current status of PISNe today is rather grim. Recent nu- 
meric al simulations indicate that super-massive stars are hard to 
form jHosokawa et alj|201 ih . If they did form, we do not currently 
expect them to exp lode as Type Ic PISNe unless at very low Z 
jUanger et al.ll2007n . If they did form and exploded at Zq/S as a 
Type Ic SN, we do not expect them to resemble SN 2007bi at all, 
as demonstrated here, but rather display red colors and lines with a 
modest width at all times. 

A more attractive scenario for producing blue weakly- 
blanketed broad-lined superluminous SNe is by the magnetar 
scenario, or any similar setup in which a large energy is fed 
into th e ejecta once it has expanded to a SN size jPessart et al.l 
l2012d) . Two effects suggest that this scenario is more amenable 
to explain these events. First, the formation of a magnetar in- 
volves the fast rotation of the progenitor Fe core prior to col- 
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lapse ( iHirschi et al1l2004bl : Iwooslev & Hege3l2006l : [Oeorgv et all 
l2009l) . which can be produced in a wide range of massive stars at 
LMC-like metalhcities or less. We expect a diversity both in the 
magnetar properties and in the associated SN ejecta. Importantly, 
a large magnetar radiation is function of the magnetar properties 
and can therefore occur in combination with low- or moderate- 
mass ejecta. In contrast, a large ^''Ni mass in a PISN requires 
by essence a larger mass progenitor The snow-plow effect from 
magnetar energy injection can thus easily alter the dynamics of 
the inner e jecta and lead to t he formation of broad lines at neb- 
ular ti mes jKasen & Bildsteni r2010: Woosley 2010: lOessart et al.l 
l2012cl) . Mistaking such a magnetar-powered ejecta for a PISN (in 
homologous expansion and showing a ^''Ni-powered LC) would 
lead to a poor estimate of the ejecta mass and kinetic energy. Such 
inconsistencies have arisen for SN2005bf a nd SN2007bi if one 
assumes a super- massive ^^ Ni-po wered SN jpolatelli et id]|2006l ; 



-Yam et al 



20091: iKasen et al. 201 ll) rather than a magnetar- 



powered SN jMaeda et al..,2007. : .Kasen & Bildste"rill2010h . What 



we demonstrate in this work is that distinguishing the two models 
should be done through detailed spectroscopic modeling. 

Our simulations give the fundamental properties we expect for 
PISNe in the early Universe and stemming from non-rotating pro- 
genitor stars. Proposed PISN candidates have all been observed in 
the nearby Universe, at metalhcities close to the solar value. At 
smaller metalhcities, the evolution of the progenitor massive star 
will be altered, in ways that depe nd considerably on the mass loss 
rate. Exploratory calculations bv lLanger et alj | |2007|) suggest that 
H-rich PISNe are possible at metalhcities as high as Zq/3, but that 
H-deficient (Type Ic) PISNe may be limited to lower-metallicity 
environments. Besides the impact on the progenitor properties (and 
its propensity to lead to core collapse or pair production), variations 
in the metallicity will tend to modulate the effects of line blan- 
keting. Based on our results, we anticipate this will play a role at 
early times (plateau phase in SNe II-P, pre-peak phase in SNe II- 
pec and Ib/c), before the photosphere recedes to the inner ejecta 
rich in explosively-synthesized IGEs. 

A potential caveat of our simulations, to be remedied in the fu- 
ture, is the neglect of stellar rotation for the pre-SN evolution. The 
PISN cha nnel for rotating massive Pop IIII star s has been recentl y 
studied bv lChatzopoulos & Wheelerl ( 1201 2h and lYoon etakl ( l2012h . 
By increasing the stellar core mass, fast rotation can lower the mini- 
mum main-sequence mass needed to enc ounter the pair-production 
instability. Neglecting mass loss entirelv- Tchatzopoulos & Wheeled 
( I2OI2I) suggest a lower mass limit of ~ 65 M© ; allowing for 
radiatio n-driven mass los s from optically-thick lines of CNO el- 
ements, |%on|eral] l|20l3) give a lower main-sequence mass limit 
of ~9OM0, corresponding to a final mass of ~7OM0. Rotation 
can thus permit to produce lower mass PISN ejecta, about 30% less 
massive than our HelOO model. Lowering the ejecta mass reduces 
the rise time to light-curve peak, although this shift will depend 
on the explosion energy and '^^m mass produced. This will di- 
versify the PISN light-curve and spectral properties, introducing a 
range o f peak luminosities and spectral properties. In lOessart et al.l 
( l2012d) . we presented results for PISNe from He cores between 
100 and 125 Mq (evolved without mass loss). These covered a huge 
range of ^^Ni mass between 5.0 and 32.4 Mq, ejecta kinetic en- 
ergy between 37.6 and 74.2 B. More energetic models produced 
higher peak luminosities with br oader spectral lines at any given 
epoch. However, as emphasized in lOessart et alj l l2012cl) . the ejecta 
richer in ^®Ni have larger temperatures but suffer stronger blan- 
keting from IGEs so that the range of colors is limited and the re- 
sulting PISN spectra are rather red. To conclude, the variations of 



a few tens of percent in ejecta mass permitted by the inclusion of 
progenitor rotation will likely have only a modest impact on PISN 
radiative properties. 

Ongoing and future deep systematic surveys of the sky may 
eventually detect unambiguously the explosion of the first stars. 
Given their longer high-brightness light curve and bluer colours, 
the explosion of RSG stars should be prime targets. Their extended 
and massive hydrogen-envelope would also allow the inference of 
the environmental metallicity in which they form. Furthermore, 
such RSG-star explosions could be used to determine distances 
out to large redshifts us i ng the Expanding Photosphere Method 
jKirshner & KwanI 1 1974 lEastman et all 1 19961 : iDessart & Hillied 
l2005ah . 
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APPENDIX A: EJECTA PROPERTIES 

In this section, we present more quantitatively the photospheric and 
ejecta properties. We give tabulated values (Tables lATl - 4A4t corre- 
sponding to Fig. 5 discussed in Section 4. Specifically, we give the 
time evolution for the radius, velocity, temperature, and overlying 
mass at the electron-scattering photosphere, together with the base 
ejecta electron-scattering optical depth. 

Concerning the latter, we note that ejecta recombination in- 
duces a decline of the optical depth that is steeper than the 1/t^ 
scaling corresponding to homologous expansion. Nonetheless, the 
electron-scattering optical depth remains > 1 for about 600 d in all 
four simulations. This is a rather unique property of PISN explo- 
sions, inherent to the huge progenitor/ejecta masses involved. This 
property may be misleading because by the time the ejecta base is 
optically thin, the outer ejecta has a much lower density and may 
thus be in a nebular regime. This hybrid configuration is reflected 
spectroscopically by the concomitant presence of P-Cygni profiles 
(associated with transitions that remain optically thick for years) 
and forbidden emission lines (such as the [Ca II] 7300 A doublet). 
There is thus continued interaction between the radiation and the 
gas, but the depletion of the radiation field at short wavelength (in 
the bound-free continua of H/He/CNO/IMEs) is very weak and thus 
unable to cause much photo-ionization. At late times, 7-rays travel 
some distance before being absorbed, inducing non-local energy 
deposition. The corresponding thermal energy, and the associated 
non-thermal excitation and ionization, may cause the ionization to 
rise again. 

We note that the huge energy release from decay (and the 
diffusion of that heat) prevent the formation of a sharp drop in 
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Table Al. Evolution of important quantities at the electron-scattering pho- 
tosphere for model B 190. By definition, these are limited to so-called photo- 
spheric epochs. They are also primarily illustrative since electron scattering 
becomes irrelevant as metal-line opacity increase and eventually dominate 
at late times. 

B190 



Age 




Kphot 




Tphot 


AMphot 


[d] 




[10^^ cm] 


[km s"^] 


[K] 


[Mq] 


16.7 


8456.35 


3.06 


21224 


5589 


0.313 


18.3 


6805.03 


3.18 


20058 


5318 


0.412 


20.2 


5436.70 


3.25 


18627 


5101 


0.563 


22.2 


4317.84 


3.23 


16870 


4863 


0.801 


24.4 


3436.37 


3.14 


14874 


4680 


1.152 


26.9 


2754.18 


2.97 


12821 


4549 


1.650 


29.5 


2226.68 


2.80 


10954 


4456 


2.485 


32.5 


1810.37 


2.66 


9462 


4411 


7.520 


35.7 


1471.41 


2.51 


8143 


4446 


13.200 


39.3 


1180.57 


2.34 


6876 


4453 


18.964 


43.2 


940.03 


2.03 


5432 


4308 


28.061 


47.6 


750.98 


1.78 


4325 


4225 


44.386 


52.3 


605.75 


1.86 


4104 


4250 


50.383 


57.6 


493.08 


1.98 


3978 


4376 


54.026 


63.3 


403.21 


2.12 


3878 


4557 


56.901 


69.7 


329.64 


2.28 


3795 


4752 


59.300 


76.6 


268.23 


2.47 


3727 


4946 


61.280 


84.3 


219.36 


2.68 


3675 


5128 


62.839 


92.7 


179.38 


2.91 


3637 


5294 


63.989 


102.0 


146.02 


3.18 


3607 


5447 


64.867 


112.0 


118.34 


3.46 


3579 


5583 


65.699 


122.0 


96.86 


3.75 


3554 


5688 


66.461 


132.0 


80.21 


4.02 


3522 


5769 


67.447 


142.0 


66.94 


4.27 


3479 


5830 


68.769 


152.0 


56.06 


4.50 


3424 


5871 


70.493 


162.0 


46.49 


4.59 


3276 


5855 


75.231 


172.0 


38.57 


4.70 


3163 


5866 


78.905 


182.0 


31.97 


4.81 


3057 


5869 


82.450 


192.0 


26.73 


4.89 


2949 


5850 


86.025 


211.0 


19.27 


5.01 


2750 


5735 


92.428 


232.0 


13.32 


5.18 


2584 


5432 


97.182 


255.0 


9.19 


5.43 


2465 


4938 


100.306 


280.0 


6.41 


5.71 


2359 


4552 


102.933 


308.0 


4.52 


5.98 


2247 


4277 


105.600 


338.0 


3.43 


6.24 


2135 


4029 


108.257 


371.0 


2.82 


6.47 


2017 


3677 


110.939 


408.0 


2.31 


6.66 


1890 


3225 


113.995 


448.8 


1.88 


6.78 


1748 


3037 


117.161 


493.7 


1.51 


6.78 


1590 


2934 


120.450 


543.0 


1.19 


6.84 


1458 


2835 


123.780 


597.0 


0.96 


6.80 


1318 


2780 


127.295 


657.0 


0.78 


5.87 


1034 


2891 


130.894 



optical depth at any time, something that is most dramatically 
seen i n "standard" SNe II-P sim ulations at the end of the plateau 
phase jPessart & Hillieill201 lij) . In that case, the transition to this 
nebular regime often coincides with a sharp rise in polarization 
( [Leonard et alllOOdlPessart & Hillieill201 Ibl) . 

We also show the ionization state of the gas in Fig. lAll as it in- 
dicates which species contribute opacity at depth (thus controling 
the diffusion of heat) and at the photosphere (affecting the spec- 
trum formation). We restrict the illustration to models R190, B210, 
and HelOOionI, and the dominant species H and/or He, O, and Fe. 
A generic property of these ejecta is their low ionization. This re- 
sult may seem counter intuitive given the large amount of ^''Ni but 



Table A2. Same as Table lAl] but now showing the photospheric properties 
for model R190. Note that in the radiative-transfer model R190 the total 
ejecta mass is 7% more massive than the original non-homologously ex- 
panding Id input at 30 d (see Footnote 2 for details). 
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5316 


4466 


33.147 


240.8 


11.69 


9.69 


4655 


4139 


38.840 


264.0 


7.30 


6.96 


3052 


4610 


130.108 


290.0 


4.70 


6.44 


2569 


4416 


144.765 


319.0 


3.31 


6.49 


2353 


4131 


150.075 


350.9 


2.67 


6.66 


2196 


3734 


153.387 


386.0 


2.18 


6.86 


2057 


3395 


156.586 


424.6 


1.78 


7.07 


1926 


3097 


159.586 


467.0 


1.44 


7.26 


1798 


2956 


163.258 


513.7 


1.15 


7.38 


1661 


2891 


166.955 


565.0 


0.92 


7.24 


1483 


2897 


170.969 


621.5 


0.75 


5.68 


1058 


3060 


175.209 



it stems from the large ejecta mass which increases the diffusion 
time, making the SN evolve on very long time scales and turning 
transparent after few hundred days when the decay energy rate is 
small. This large mass also means this decay energy is to be shared, 
so that the energy released per unit mass is in fact quite compara- 
ble to what obtains in standard core-collapse SNe. Combined with 
the huge explosion energies, the small initial radius implies a very 
large expansion, associated with strong cooling from PdV work. 

These curves reproduce the general morphology of ionization 
profiles in "standard" SN ejecta (e.g., Dessart & Hillier 2010.). The 
fast-expanding low-density outer regions maintain a high ioniza- 
tion - this ionization freeze-ou t is in a large part a time-dependent 
effect jPessart & HilhedbOOSh . At low velocity, radioactive decay 
contributes significant heating, weakly affected by radiative cool- 
ing due to the larger optical depth of the inner ejecta layers. The 
photosphere (marked as a filled circle in the figure) is by essence 
a tracer of the region bridging thick and thin conditions, and thus 
delimits the higher ionization ejecta at depth from the recombined 
conditions immediately above it. The minimum ionization is found 
at the photosphere up to the peak of the light curve because it suf- 
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Figure Al. Evolution of the ionization state of selected species versus time (color coding) and velocity (equivalent to the depth in the ejecta) for models 
R190NL (top row; we show H, O, and Fe), B210NL (middle row; He, O, and Fe), and HelOONL (bottom row; He, O, and Fe). An ionization level n, where n 
is a positive real number, for a species X means that this species X is found at that location primarily in its n"^-time ionized form (i.e., the ionization level is 
equal to X'+ ). Along each curve, a dot represents the photospheric location at each epoch. 



fers str ong radiative cooling, fast expansion cooling, and weaker 
heatingEj 

Finally, we illustrate in Fig. lA2l the temperature evolution of 
the B210 model ejecta, and in particular the presence of a rever- 
sal in the trajectory of the photosphere in velocity/mass space. 
This reversal is also seen in the ionization of oxygen, for example 
(Fig. lAU . This occurs here because of the strong heating from de- 
cay. Early on, the photosphere recedes in mass until the heat wave 
from greater depth halts its recession, and in the present (but rare) 
situation, reverses that recession to outward migration. We obtained 
such an effect in SNe Ib /Ic models in which no macroscopic mix- 
ing of ^''Ni was applied jPessart et al.ll201 ih . Observationally, this 
may occur if a strong heat source is turned on, as in a magnetar. 
SN 2005bf may well have been influenced by such a delayed heat- 



We note that species tend to be over-ionized when subdominant so these 
plots have to be interpreted with the chemical stratification in mind (Fig. 2). 
We also note a peculiar dip in ionization at 2000-4000 km s^^ in model 
R190, which is associated with the O-rich shell: This shell is highly-bound 
but ^^Ni deficient so that at early times when diffusion is inhibited the cor- 
responding layers evolve adiabatically and cool tremendously by expanding 
from their small original radius. 



ing, given the observed double-peaked light curve, the reversal of 
the velocity at maximum absorpt ion in He 1 5875 A , and the blue 
colors of the spectra at late times jMaeda et alJl2007l ; lDessart et al.l 
l2012cl) . 



APPENDIX B: ATOMIC DATA SOURCES AND MODEL 
ATOMS 

The model atoms adopted for all simulation s in this work are identi- 
cal to that used in iDessart & Hillieil l l201 lij) . The sources of atomic 
data are varied, and in many cases multiple data sets for a given ion 
are available. In some cases these multiple data sets represent an 
evolution in data quality and/or quantity, while in other cases they 
represent different sources and/or computational methods. Com- 
parisons of models calculated with different data sets and atomic 
models potentially provide insights into the sensitivity of our re- 
sults to the adopted model atoms and hydrodynamical inputs (al- 
though such calculations have yet to be undertaken for SNe). 

Oscillator strengths for CNO elements were originally taken 
from iNussbaumer & StorevI ( Il983l , Il984l) . These authors also pro- 
vide transition probabilities to states in the ion continuum. The 
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Table A3. Same as Table lAl] but now showing the photospheric properties 
for model B210. 



Table A4. Same as TableEH but now showing the photospheric properties 
for model He 100. 



B210 



HelOO 



Age 


^basc.cs 


^■phot 


^phot 


^phot 




[d] 




[10^^ cm 


] [kms-1] 


[K] 


[Mq] 


16.8 


3806.45 


3.08 


21214 


5453 


0.491 


18.5 


3090.62 


3.12 


19555 


5155 


0.605 


20.3 


2496.34 


3.07 


17476 


4934 


0.786 


22.4 


2009.13 


2.98 


15414 


4741 


1.362 


24.6 


1631.11 


2.93 


13793 


4645 


2.656 


27.1 


1329.80 


2.90 


12380 


4601 


5.131 


29.8 


1085.34 


2.84 


11041 


4582 


7.947 


32.8 


876.82 


2.79 


9849 


4565 


11.126 


36.1 


707.42 


2.79 


8951 


4550 


14.555 


39.7 


572.21 


2.53 


7373 


4704 


26.563 


43.6 


461.90 


2.65 


7031 


5102 


35.364 


48.0 


373.02 


2.87 


6922 


5421 


39.163 


52.8 


302.94 


3.15 


6894 


5699 


40.280 


58.1 


246.63 


3.47 


6920 


5935 


39.177 


63.9 


199.83 


3.87 


7002 


6117 


36.066 


70.3 


160.92 


4.34 


7152 


6245 


31.235 


77.3 


129.92 


4.94 


7393 


6360 


25.803 


85.0 


105.40 


5.66 


7695 


6526 


21.755 


93.6 


85.15 


6.41 


7929 


6687 


19.797 


102.0 


69.46 


7.08 


8038 


6776 


19.035 


112.0 


55.20 


7.78 


8038 


6823 


19.055 


122.0 


44.78 


8.36 


7933 


6820 


19.792 


132.0 


36.76 


8.81 


7728 


6715 


21.504 


142.0 


30.26 


9.20 


7495 


6589 


24.283 


156.0 


22.26 


9.68 


7181 


6387 


30.617 


171.0 


17.00 


10.15 


6873 


6242 


41.190 


188.0 


13.15 


10.49 


6455 


6264 


58.370 


207.0 


9.81 


10.33 


5777 


6502 


80.717 


227.0 


6.83 


9.51 


4849 


6785 


100.010 


249.0 


4.90 


9.23 


4288 


6405 


110.797 


273.9 


3.79 


9.40 


3971 


5895 


117.728 


301.3 


2.81 


9.47 


3639 


5578 


124.344 


331.4 


1.77 


8.98 


3136 


5502 


130.862 


364.5 


1.42 


7.92 


2515 


5454 


136.697 


401.0 


1.16 


6.78 


1956 


5294 


141.143 


441.1 


0.95 


5.46 


1432 


5099 


144.326 


485.2 


0.78 


3.69 


881 


4880 


146.294 


533.7 


0.63 


2.00 


433 


4621 


146.849 


587.1 


0.51 


2.20 


433 


4320 


146.833 



largest source of oscillator data is from iKurucj l l2009l |2010|) ; 
its principal advantage over many other sources (e.g., Opacity 
Project) is that LS coupling is not assumed. More recently, non- 
LS osci llator strengths have become available through the Iron 
Project jHummer et al.|[l993h . an d work done by the atomic-data 
group at Ohio State University (iNahatI lOlOh. Ot her important 
source s o f radiat i ve dat a for Fe include iBecker & B utler ( 199^ 
Il995al lbl). iNahaJ i ll995l) . Atomic data from the opacity project 
comes from TOPBASE l lCunto et al.l Il993l) . Energy levels have 
generally been obtained from National Institute of Standards and 
Technology. Collisional data is sparse, particularly for states far 
from the ground state. The principal source for collisional data 
a mong low l ying s tates for a variety of species is the tabulatio n 
by Mendozaj (Il983h: ot her sources includ e lBerrington et al] 1 19851). 
iLennon et all ill 9851) . [Lennon & Burk3 ( Il994b IShine & Linskvl 
( ll974l) . lTava1 l ll997al M. Zhang & Pradhan ( Il995al : ll995bl : ll997l) 
Photoionization data is taken from the Opacity Project jSeatonl 



Age 


^basc.cs 


^ptiot 


Vphot 


^phot 


AMphot 


[d] 




[lO'^^ cm] 


[km s^i] 


[K] 


[Mq] 


11.6 


6433.32 


0.90 


8977 


5183 


8.295 


12.7 


5215.37 


0.96 


8746 


4992 


9.262 


14.0 


4236.04 


1.03 


8535 


4926 


10.345 


15.4 


3449.59 


1.10 


8284 


4861 


11.927 


16.9 


2816.45 


1.17 


7979 


4791 


14.178 


18.6 


2300.95 


1.23 


7650 


4742 


16.979 


20.5 


1875.86 


1.29 


7319 


4699 


20.203 


22.5 


1524.64 


1.37 


7046 


4598 


23.147 


24.8 


1237.52 


1.47 


6855 


4435 


25.342 


27.2 


1005.06 


1.58 


6726 


4275 


26.909 


30.0 


818.25 


1.71 


6604 


4207 


28.425 


33.0 


666.34 


1.85 


6475 


4237 


30.085 


36.3 


541.71 


1.99 


6351 


4327 


31.746 


39.9 


438.67 


2.14 


6218 


4465 


33.583 


43.9 


355.25 


2.31 


6084 


4621 


35.513 


48.3 


288.98 


2.48 


5953 


4800 


37.467 


53.1 


235.24 


2.68 


5828 


4984 


39.364 


58.4 


190.65 


2.90 


5746 


5166 


40.677 


64.3 


153.97 


3.16 


5683 


5333 


41.662 


70.7 


124.85 


3.47 


5685 


5470 


41.650 


77.8 


101.46 


3.82 


5684 


5597 


41.661 


85.6 


81.86 


4.22 


5706 


5704 


41.316 


94.1 


65.21 


4.64 


5700 


5804 


41.418 


103.5 


51.21 


4.95 


5532 


5876 


44.158 


113.0 


41.30 


5.28 


5409 


5954 


46.228 


124.3 


32.20 


5.62 


5235 


6022 


49.186 


136.7 


24.34 


5.91 


5002 


6075 


53.297 


150.4 


17.93 


6.10 


4691 


6150 


58.696 


165.0 


13.51 


6.15 


4314 


6189 


64.763 


181.0 


10.05 


6.20 


3962 


6055 


69.968 


199.0 


7.04 


6.35 


3691 


5704 


73.820 


219.0 


5.10 


6.56 


3464 


5286 


76.981 


240.9 


3.80 


6.74 


3236 


4953 


80.037 


265.0 


2.79 


6.86 


2995 


4750 


83.083 


291.5 


2.05 


6.90 


2741 


4641 


86.162 


320.6 


1.68 


6.89 


2488 


4575 


89.257 


352.7 


1.38 


6.80 


2230 


4518 


92.417 


388.0 


1.13 


6.44 


1919 


4466 


95.531 


426.8 


0.93 


5.47 


1482 


4435 


98.325 


469.5 


0.76 


3.68 


906 


4398 


100.230 


516.5 


0.62 


2.03 


455 


4240 


100.789 


568.1 


0.50 


2.23 


455 


3970 


100.789 


624.9 


0.40 


2.46 


455 


3701 


100.789 



I987I: ICuntoetalJ|l993l) . the Iron Project jHtimmer et al.l [19931 : 
Nahar & PradhaiJ (l996h . and iNahar & PradhanI ( ll993l) rUnfortu- 



nately Ni and Co photoionization data is generally unavailable, 
and we have utilized crude approximations . Charge exchange 
cross- sections are from the tabulation by iKingdon & Ferland 
199^ ). Atomic data for C IV was obtained from iLeibowitz 



19721) : [Peach et alj lll988l). and for the carbon isoelectronic se- 
quence from Luo & Pradhan (19891). Collision strengths for Aril 
are from |Tayal & Henry ( 1996J. The LS Nel photoioni zation 
cross-sections were modified according to ISeatonI l ll998l) . The 
same procedure was applied to using An mixing coefficients 
computed at http://aphysics2.lanl.gov/tempweb/lanl Additional 
data for Ne I was obtained from the MCHF/MCDHF web site: 
|http://nlte.mst.gov/MCHF, 
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Table Bl. Summary of the model atoms used in our PISN calculations with 
CMFGEN (H is not included in model HelOO). Nf refers to the number of 
full levels, Ns to the number of super levels, and Ntrans to the correspond- 
ing number of bound-bound transitions. The last column refers to the upper 
level for each ion treated. At late times, we exclude the high ionization 
stages and simultaneously increase the number of levels for Fe I and Fe II. 
After a few hundred days, we also split the lower 5, 10, or 15 super-levels 
to account explicitly for more processes between low-lying states of IGEs. 
nw refers to a state with principal quantum number n (all I states com- 
bined into a single state), and spin 2. Similarly, 8z refers to the n = 8 
state with high 1 states (usually I — 4 and above) combined and spin 1 . 
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Figure A2. Evolution of the ejecta temperature for model B210 from 25 
until ~ 600 d after explosion. The dots refer to the location of the electron- 
scattering photosphere and highlight the non-monotonic behavior of the 
photospheric location within the ejecta. Indeed, the strong heat wave pow- 
ered by the decay energy from 21.3 Mq of ^®Ni causes the photosphere to 
migrate outwards in mass (or velocity) between 50 and lOOd after explo- 
sion. This behavior, unique to model B210, is seen in both mass and velocity 
at the photosphere in Fig. 5. Note that despite the huge ^^Ni mass, the tem- 
perature at a few hundred days after explosion is < 5000 K throughout the 
ejecta. 



APPENDIX C: DEPENDENCY ON MODEL ATOMS 

The main impact of increasing the size of the model atoms in our 
simulations, in particular for metals, is to enhance the magnitude 
of line blanketing, which tends to make the SED redder. This effect 
is generally weaker than obtained here through variations in com- 
position between models or versus time (i.e., H vs. IMEs, IMEs vs. 
IGEs etc.), although it noticeably alters the colors after the light- 
curve peak. 

For the early-time simulations of models R190, B190, and 
B210, the photosphere is located in essentially a pure Hydro- 
gen/Helium plasma. Line blanketing effects are negligible and there 
is no concern with the opacity of IMEs and IGEs and the complete- 
ness of our model atoms at such early times. 

However, as the photosphere recedes to deeper layers in 
the ejecta where IMEs and IGEs are abundant the complete- 
ness of model atoms becomes a central concern for the reliabil- 
ity of our predictions. After much experimentation in past simula- 
tions toessart & Hilliei|[201ol l2011a| - iLi et"ai]|2012l : iDessart et al] 
l2012dir we have converged to an adequate assignment for the num- 
ber of full and super le vels to include (for a discussion on our 
super-level approach, see lHillier& Millet! 1998h . Recently, we em- 
phasize d the critical n eed of including Fe I in simulations of Type 
II SNe jLi et al.ll2012l) . while earlier on we found that Sc and Ti 
play a critical role in optical spectra of SNe II-P jDessart & Hillied 
l2011al) . 

In ejecta dominated by metals, the situation is more compli- 
cated than in type II SNe. To reach a reliable radiative-transfer re- 
sult, we first find that more levels need to be included. Secondly, 
while we typically include in the non-LTE solver all the metal line 
transitions with a statistical weight grea ter than 0.002, we find that 
we now need to go down to lO^''!^ This means treating many 



This cut only applies to elements beyond Ne in the periodic table, does 
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more lines in total, which increases the computation time. Finally, 
in the course of this study started with models R190, B190, B210, 
and HelOO, we found that the ejecta were rather cold at the nebular 
phase, and for the compact progenitors, these conditions were also 
cold prior to re-brightening at early times. Because of the domi- 
nance of IMEs in the corresponding emitting layers, we have run 
a companion model, named HelOOionI, identical to HelOO (which 
has the same model atoms as R190 etc), but now including the neu- 
tral states for species Mg, Si, S, and Ca. We compare these two 
sequences at selected post-explosion times in Fig. lCll 

The bolometric light curve is largely insensitive to the addi- 
tion of these neutral atoms. We understand this from the fact that the 
light curve is primarily conditioned by the diffusion of heat through 
the optically-thick layers of the ejecta. These shells tend to be ion- 
ized and hot, and thus dominated by more ionized species. The only 
slight change occurs at very early times, prior to re-brightening, 
when the photosphere is very cold (i.e., ~4000K; Fig. 5). At such 
times, the blanketing from neutral species is actually huge, because 
of the IME-rich photosphere, but also because the IGEs are under 
abundant compared to their solar-metallicity value. Hence, these 
neutral species act as an overlying blanketing layer, modifying the 
blocking power of the last-scattering layer. The blanketing is in part 
due to lines, especially in the optical and the near-IR, but more im- 
portantly stems from bound-free opacity from the ground and ex- 
cited states (note in particular the photo-ionization edge of Mg I at 
3757 A) of the corresponding ions (see also the continuum curve in 
Fig. 14). As the photosphere heats up on the rise to peak, species 
become once-ionized and the difference between models HelOO 
and HelOOionI weakens, essentially limited to the presence of the 
strong Mg 1 5 173 A line. 

However, as the SN light curve passes its peak, the color 
reddens again and the photosphere cools down to 4000 K within 
100-200 d. The SN becomes nebular and blanketing from neutral 
species turns back on. At nebular times, additional blanketing red- 
dens the SED below 5500 A, but also impacts the long wavelength 
range. This stems indirectly from the impact on the blanketed ra- 
diation field. It also arises directly from additional optically-thick 
lines from Mg I, Si I or S I. 

As discussed above, metal line blanketing is irrelevant up to 
the peak because the composition really becomes IGE dominated 
at the photosphere only then. So, to expedite our very time con- 
suming simulations, we have used modest-sized iron model atoms 
up to the peak (Fe I [136,44] and Fe II [1 15,50]), improved it early 
after the peak (Eel [136,44] and Fell [827,275]), and finally used 
huge model atoms for both at later times (Eel [1142,413] and 
Fell [827,275]) — the numbers in square brackets represent the 
number of full and super levels. This choice of model atoms yields 
converged results in our SNe la simulations l lDessart et al.ll2012ah 
and we expect the same level of accuracy to be reached in those 
similarly metal-rich PISN ejecta at late times. 



not apply to the lowest n levels (n is typically 9), and we only cut a transi- 
tion when there is at least m (m is typically 9) stronger downward transi- 
tions from the level. Thus, this procedures does not cut important transitions 
to ground levels, and forbidden and semi-forbidden transitions among low- 
lying states. 



APPENDIX D: LOG OF MAGNITUDES FOR EACH PISN 
MODEL 

To complement the discussion on the colors of our PISN simula- 
tions, we present in Tables lDTl4D4l the magnitudes and luminosities 
(bolometric and UVOIR, i.e., integrated flux from the blue edge of 
the U band to the red edge of the I band). From these, one can easily 
infer various bolometric corrections and compare to observations. 
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Figure CI. Illustration of the changes to the bolometric light cui-ve (top left) and the spectra (subsequent panels, showing snapshots at 30, 50, 100, 200, 400, 
and 620 d after explosion) for the ejecta model HelOOionlNL, which treats the neutral species Mg I, Ca I, Si I, and S I, with model HelOONL, which ignores 
them. Under favorable ionization conditions, these additional neutral states can increase the magnitude of Hne blanketing and bound-free opacity. Notice at 
nebular times the impact on the strength of the [Call] 7300 A doublet (HelOONL), reduced by the overlap with Cal Unes at 7148 and 7326A (HelOOionNL). 
Opacity effects are associated with line blanketing and photo-ionization cross sections, and are strongly modulated by the ionization state at the photosphere. 
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Table Dl. Bolometiic and UVOIR luminosities as well as synthetic magnitudes for PISN model B 190NL. The UVOIR luminosity coiTesponds to the integrated 
flux over the wavelength interval 3000-9200 A, i.e., from the blue edge of the U band to the red edge of the I band. 



Age Lboi iuvoiR Mbol Mu Mb Mv Mr Mi Mj Mh Mk 

(d) (ergs-i) (ergs-l) (mag) (mag) (mag) (mag) (mag) (mag) (mag) (mag) (mag) 

16.67 3.4929e+42 2.2666e+42 -17.651 -16.568 -16.938 -17.578 -18.068 -18.129 -18.560 -18.677 -18.926 

18.34 3.0239e-i-42 1.9194e-i-42 -17.495 -16.028 -16.691 -17.401 -17.951 -18.011 -18.506 -18.618 -18.900 

20.17 2.6104e-i-42 1.6095e-i-42 -17.335 -15.449 -16.419 -17.207 -17.814 -17.876 -18.429 -18.534 -18.849 

22.19 2.1978e-i-42 1.3107e-i-42 -17.148 -14.824 -16.111 -16.960 -17.638 -17.701 -18.310 -18.408 -18.755 

24.41 1.8688e-i-42 1.0867e-i-42 -16.972 -14.222 -15.816 -16.764 -17.464 -17.540 -18.176 -18.274 -18.633 

26.85 1.6860e-i-42 9.7437e-i-41 -16.860 -13.851 -15.654 -16.661 -17.344 -17.458 -18.068 -18.187 -18.528 

29.54 1.5034e-i-42 8.6590e-i-41 -16.736 -13.739 -15.543 -16.501 -17.210 -17.338 -17.935 -18.081 -18.403 

32.49 1.3165e-i-42 7.5756e-i-41 -16.592 -13.666 -15.430 -16.344 -17.058 -17.181 -17.784 -17.944 -18.254 

35.74 1.1436e-i-42 6.6064e-i-41 -16.439 -13.570 -15.324 -16.189 -16.901 -17.018 -17.624 -17.788 -18.093 

39.32 9.6780e-t-41 5.7329e-H41 -16.258 -13.577 -15.206 -16.009 -16.739 -16.842 -17.428 -17.480 -17.826 
43.25 7.5661e-t-41 4.7465e+41 -15.990 -13.855 -14.965 -15.849 -16.531 -16.629 -17.056 -17.097 -17.410 
47.58 6.6828e+41 4.3700e-H41 -15.856 -14.454 -14.928 -15.837 -16.355 -16.518 -16.850 -16.976 -17.177 
52.34 7.2292e-i-41 4.8332e-i-41 -15.941 -14.920 -15.111 -15.981 -16.380 -16.608 -16.907 -17.071 -17.229 
57.57 9.0472e-i-41 6.1748e-i-41 -16.184 -15.382 -15.452 -16.244 -16.584 -16.840 -17.113 -17.286 -17.426 

63.33 1.2072e-i-42 8.4140e-i-41 -16.498 -15.860 -15.854 -16.561 -16.870 -17.130 -17.373 -17.551 -17.682 
69.66 1.6449e-i-42 1.1691e-i-42 -16.833 -16.338 -16.261 -16.897 -17.185 -17.436 -17.652 -17.830 -17.958 
76.63 2.2447e-i-42 1.6236e-i-42 -17.171 -16.804 -16.656 -17.231 -17.504 -17.740 -17.931 -18.108 -18.235 
84.29 3.0375e-i-42 2.2306e-i-42 -17.499 -17.248 -17.031 -17.554 -17.815 -18.033 -18.203 -18.375 -18.503 
92.72 4.0579e-i-42 3.0184e-i-42 -17.814 -17.664 -17.383 -17.859 -18.110 -18.314 -18.464 -18.630 -18.758 
102.00 5.3415e-i-42 4.0149e-i-42 -18.112 -18.050 -17.712 -18.146 -18.389 -18.581 -18.715 -18.870 -18.997 
112.00 6.8734e-i-42 5.2077e-i-42 -18.386 -18.394 -18.009 -18.408 -18.645 -18.829 -18.947 -19.090 -19.215 
122.00 8.5088e-i-42 6.4832e-i-42 -18.618 -18.679 -18.258 -18.630 -18.861 -19.042 -19.147 -19.278 -19.398 
132.00 1.0182e-i-43 7.7896e-i-42 -18.813 -18.911 -18.466 -18.817 -19.045 -19.223 -19.319 -19.437 -19.551 
142.00 1.1816e-i-43 9.0680e-i-42 -18.974 -19.098 -18.637 -18.973 -19.199 -19.378 -19.464 -19.572 -19.677 
152.00 1.3333e-i-43 1.0257e-i-43 -19.105 -19.241 -18.776 -19.103 -19.328 -19.508 -19.585 -19.685 -19.780 
162.00 1.4574e+43 1.1338e-i-43 -19.202 -19.264 -18.910 -19.280 -19.476 -19.642 -19.672 -19.772 -19.853 
172.00 1.5608e-i-43 1.2171e-i-43 -19.277 -19.290 -18.988 -19.374 -19.576 -19.740 -19.748 -19.847 -19.916 
182.00 1.6339e-i-43 1.2731e-i-43 -19.326 -19.271 -19.029 -19.442 -19.654 -19.824 -19.808 -19.908 -19.968 
192.00 1.6723e-i-43 1.2992e-i-43 -19.351 -19.206 -19.033 -19.485 -19.712 -19.891 -19.849 -19.952 -20.007 
211.00 1.6321e-i-43 1.2504e-i-43 -19.325 -18.946 -18.920 -19.485 -19.755 -19.956 -19.873 -19.999 -20.051 
232.00 1.4259e-i-43 1.0593e-i-43 -19.178 -18.465 -18.617 -19.341 -19.668 -19.902 -19.822 -19.998 -20.049 
255.00 1.0714e-l-43 7.5083e-l-42 -18.868 -17.487 -18.050 -18.991 -19.403 -19.687 -19.684 -19.894 -19.951 
280.00 7.4220e-i-42 4.9189e-i-42 -18.469 -16.488 -17.413 -18.480 -19.008 -19.372 -19.417 -19.659 -19.726 
308.00 5.1085e-i-42 3.2255e-i-42 -18.064 -15.741 -16.787 -17.925 -18.575 -19.050 -19.170 -19.355 -19.420 
338.00 2.9760e-i-42 1.9662e-i-42 -17.477 -15.067 -16.157 -17.303 -17.993 -18.651 -18.520 -18.585 -18.589 
371.00 2.0821e-i-42 1.1506e-i-42 -17.089 -12.730 -14.573 -16.445 -17.540 -18.301 -18.504 -18.796 -18.364 
408.00 1.5414e-i-42 7.0200e-i-41 -16.763 -11.737 -13.546 -15.679 -16.950 -17.896 -18.429 -18.752 -18.164 
448.80 1.1270e-i-42 4.3687e-i-41 -16.423 -11.322 -12.889 -15.005 -16.371 -17.447 -18.214 -18.598 -18.050 
493.70 7.7498e-i-41 3.0059e-i-41 -16.016 -11.055 -12.425 -14.518 -15.943 -17.060 -17.850 -18.163 -17.612 
543.10 4.9505e-i-41 2.1863e-i-41 -15.530 -10.741 -11.976 -14.087 -15.596 -16.734 -17.380 -17.374 -16.675 
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Table D2. Bolometric and UVOIR luminosities as well as synthetic magnitudes for PISN model R190NL. 



Age Lbol Lvvom Afbol Afy Mb Mv Mr Mj Mj Mh Mk 

(d) (ergs-i) (ergs^i) (mag) (mag) (mag) (mag) (mag) (mag) (mag) (mag) (mag) 

36.72 1.3795e+44 4.3057e+43 -21.642 -21.330 -20.380 -20.354 -20.440 -20.374 -20.329 -20.346 -20.454 

40.39 1.0017e-H44 3.6363e-H43 -21.295 -21.066 -20.193 -20.247 -20.348 -20.317 -20.317 -20.348 -20.469 

44.43 7.6372e-H43 3.2061e-H43 -21.000 -20.843 -20.055 -20.176 -20.297 -20.292 -20.333 -20.375 -20.507 

48.87 6.1774e-H43 2.9380e-H43 -20.770 -20.660 -19.958 -20.136 -20.279 -20.292 -20.368 -20.418 -20.558 

53.76 5.2790e+43 2.7743e-H43 -20.600 -20.509 -19.894 -20.116 -20.283 -20.308 -20.415 -20.472 -20.617 

59.14 4.7693e+43 2.6861e-H43 -20.489 -20.390 -19.857 -20.114 -20.304 -20.338 -20.472 -20.536 -20.685 

65.05 4.4163e+43 2.6210e-H43 -20.406 -20.278 -19.827 -20.112 -20.328 -20.370 -20.534 -20.606 -20.758 

71.56 4.1327e-i-43 2.5494e-H43 -20.334 -20.160 -19.790 -20.101 -20.347 -20.401 -20.599 -20.684 -20.841 

78.72 3.7587e-H43 2.4079e-H43 -20.231 -19.984 -19.716 -20.056 -20.341 -20.411 -20.654 -20.755 -20.920 

86.59 3.3955e-H43 2.2317e-H43 -20.120 -19.754 -19.616 -19.991 -20.320 -20.413 -20.702 -20.822 -20.999 

95.25 3.0894e-H43 2.0488e-H43 -20.018 -19.491 -19.502 -19.909 -20.284 -20.401 -20.738 -20.878 -21.073 

104.00 2.8994e-H43 1.9156e-H43 -19.949 -19.300 -19.412 -19.825 -20.248 -20.374 -20.749 -20.900 -21.122 

114.00 2.5779e+43 1.6759e-H43 -19.821 -18.996 -19.247 -19.662 -20.147 -20.265 -20.703 -20.855 -21.111 

124.00 2.1881e-i-43 1.3864e-H43 -19.643 -18.562 -19.006 -19.449 -19.995 -20.101 -20.616 -20.759 -21.044 

136.00 1.8206e-H43 1.1264e-H43 -19.444 -18.040 -18.722 -19.235 -19.835 -19.906 -20.491 -20.605 -20.933 

150.00 1.5858e-H43 9.8104e-H42 -19.294 -17.617 -18.512 -19.131 -19.735 -19.766 -20.364 -20.443 -20.791 

165.00 1.6326e-H43 1.0444e-H43 -19.325 -17.749 -18.571 -19.250 -19.793 -19.828 -20.334 -20.419 -20.724 

181.00 1.8889e+43 1.2586e-i-43 -19.484 -18.362 -18.813 -19.461 -19.931 -19.990 -20.388 -20.483 -20.746 

199.00 2.2135e-H43 1.5332e-H43 -19.656 -18.908 -19.058 -19.673 -20.080 -20.155 -20.478 -20.520 -20.749 

218.90 2.4176e-H43 1.7416e-H43 -19.752 -19.229 -19.179 -19.816 -20.181 -20.312 -20.509 -20.434 -20.652 

240.80 2.1726e-H43 1.6053e-H43 -19.636 -19.020 -19.031 -19.747 -20.132 -20.327 -20.306 -20.272 -20.420 

264.00 1.4704e+43 1.0458e-t-43 -19.212 -18.029 -18.347 -19.314 -19.788 -20.038 -19.960 -20.107 -20.212 

290.00 8.4332e+42 5.6610e-i-42 -18.608 -16.730 -17.494 -18.622 -19.171 -19.543 -19.522 -19.757 -19.820 

319.00 4.6253e-i-42 3.0791e-H42 -17.956 -15.695 -16.669 -17.858 -18.506 -19.064 -18.915 -19.104 -19.060 

350.90 2.8228e-i-42 1.7351e-H42 -17.420 -13.735 -15.362 -16.999 -17.978 -18.674 -18.665 -18.840 -18.545 

386.00 1.9654e+42 1.0575e-H42 -17.027 -12.129 -14.173 -16.242 -17.434 -18.289 -18.551 -18.723 -18.176 

424.60 1.3801e+42 6.3219e-t-41 -16.643 -11.844 -13.372 -15.471 -16.804 -17.823 -18.352 -18.624 -17.880 

467.00 9.9485e-H41 3.9880e-H41 -16.288 -11.553 -12.780 -14.843 -16.263 -17.359 -18.089 -18.417 -17.823 

513.70 6.4897e-H41 2.8621e-H41 -15.824 -11.402 -12.400 -14.435 -15.894 -17.012 -17.635 -17.752 -17.030 

565.00 4.1263e-H41 2.0683e-H41 -15.332 -11.179 -11.974 -13.990 -15.546 -16.680 -17.128 -16.889 -15.982 

621.50 2.5179e+41 1.4186e-H41 -14.796 -10.926 -11.484 -13.454 -15.148 -16.295 -16.542 -15.925 -14.811 

683.70 1.4502e-H41 9.0849e-H40 -14.197 -10.580 -10.909 -12.810 -14.678 -15.842 -15.859 -14.876 -13.560 

752.10 7.8231e-H40 5.5235e-H40 -13.527 -10.144 -10.213 -12.202 -14.167 -15.324 -15.121 -13.415 -11.947 

827.30 3.9231e+40 2.9709e-H40 -12.777 -9.623 -9.482 -11.362 -13.492 -14.689 -14.243 -12.345 -10.831 

910.00 1.8103e-i-40 1.4445e-H40 -11.938 -8.989 -8.676 -10.413 -12.699 -13.940 -13.244 -11.334 -9.784 

1001.00 7.6570e-i-39 6.2713e-H39 -11.003 -8.277 -7.828 -9.424 -11.785 -13.052 -12.178 -10.366 -8.765 
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Table D3. Bolometric and UVOIR luminosities as well as synthetic magnitudes for PISN model B210NL 



Age Lboi LuvoiR Mboi Mu Mb Mv Mr A/j Mj Mh Mk 

(d) (ergs~i) (ergs~i) (mag) (mag) (mag) (mag) (mag) (mag) (mag) (mag) (mag) 

16.80 2.8082e+42 1.7870e+42 -17.414 -15.893 -16.499 -17.346 -17.905 -17.961 -18.404 -18.519 -18.822 
18.48 2.4868e-i-42 1.5465e-i-42 -17.282 -15.344 -16.264 -17.196 -17.785 -17.864 -18.346 -18.464 -18.783 
20.33 2.2713e-i-42 1.3845e-i-42 -17.184 -14.879 -16.090 -17.080 -17.677 -17.792 -18.295 -18.416 -18.740 
22.36 2.1097e-i-42 1.2687e-i-42 -17.104 -14.625 -15.982 -16.974 -17.582 -17.719 -18.240 -18.365 -18.688 
24.60 1.9644e-i-42 1.1713e-i-42 -17.026 -14.490 -15.904 -16.881 -17.491 -17.639 -18.178 -18.308 -18.622 
27.10 1.8524e-i-42 1.1017e-i-42 -16.962 -14.421 -15.865 -16.813 -17.413 -17.574 -18.121 -18.259 -18.559 

29.81 1.7287e-l-42 1.0242e-l-42 -16.887 -14.350 -15.817 -16.721 -17.329 -17.495 -18.051 -18.202 -18.491 
32.79 1.7152e-i-42 1.0362e-i-42 -16.879 -14.460 -15.856 -16.762 -17.325 -17.494 -18.031 -18.141 -18.421 
36.07 2.0198e-i-42 1.2960e-i-42 -17.056 -15.322 -16.131 -17.068 -17.527 -17.684 -18.110 -18.214 -18.452 
39.68 2.6317e-i-42 1.7885e-i-42 -17.344 -16.396 -16.587 -17.421 -17.777 -17.924 -18.263 -18.408 -18.585 
43.65 3.6298e-i-42 2.5640e-i-42 -17.693 -17.200 -17.073 -17.773 -18.064 -18.213 -18.499 -18.659 -18.807 
48.01 5.0858e-i-42 3.6894e-i-42 -18.059 -17.824 -17.547 -18.121 -18.370 -18.517 -18.768 -18.931 -19.063 
52.81 7.1644e-i-42 5.2956e-i-42 -18.431 -18.370 -17.999 -18.467 -18.687 -18.825 -19.048 -19.208 -19.329 
58.09 9.9973e-i-42 7.4758e-i-42 -18.793 -18.861 -18.417 -18.797 -18.997 -19.121 -19.325 -19.475 -19.590 
63.90 1.3780e-i-43 1.0351e-i-43 -19.141 -19.307 -18.797 -19.111 -19.298 -19.407 -19.597 -19.732 -19.843 
70.29 1.8721e-i-43 1.4037e-i-43 -19.474 -19.711 -19.141 -19.408 -19.587 -19.686 -19.864 -19.980 -20.085 
77.32 2.5052e-i-43 1.8614e-i-43 -19.790 -20.067 -19.452 -19.690 -19.866 -19.954 -20.126 -20.223 -20.319 
85.05 3.3223e-i-43 2.4331e-i-43 -20.097 -20.394 -19.736 -19.962 -20.142 -20.219 -20.382 -20.460 -20.542 
93.56 4.3745e-i-43 3.1614e-i-43 -20.395 -20.719 -20.009 -20.228 -20.412 -20.480 -20.624 -20.683 -20.751 
102.00 5.5049e-i-43 3.9495e-i-43 -20.645 -20.995 -20.241 -20.454 -20.639 -20.706 -20.827 -20.868 -20.927 
112.00 6.8640e-i-43 4.9270e-i-43 -20.885 -21.263 -20.477 -20.682 -20.865 -20.934 -21.026 -21.051 -21.101 
122.00 8.1394e-i-43 5.8866e-i-43 -21.070 -21.471 -20.671 -20.868 -21.051 -21.123 -21.193 -21.200 -21.242 
132.00 9.0999e-i-43 6.7721e-i-43 -21.191 -21.626 -20.831 -21.024 -21.206 -21.278 -21.330 -21.316 -21.356 
142.00 9.9900e-i-43 7.5971e-i-43 -21.292 -21.747 -20.968 -21.154 -21.336 -21.404 -21.440 -21.403 -21.439 
156.00 1.0871e-i-44 8.5119e-i-43 -21.384 -21.850 -21.118 -21.296 -21.473 -21.532 -21.544 -21.468 -21.503 
171.00 1.1360e-i-44 9.1407e-i-43 -21.432 -21.873 -21.240 -21.418 -21.577 -21.621 -21.591 -21.476 -21.497 
188.00 1.1249e-i-44 9.2770e-i-43 -21.421 -21.757 -21.304 -21.530 -21.658 -21.674 -21.552 -21.420 -21.399 
207.00 1.0145e-i-44 8.4987e-i-43 -21.309 -21.362 -21.173 -21.563 -21.721 -21.725 -21.387 -21.329 -21.209 
227.00 8.2998e-i-43 6.8282e-i-43 -21.091 -20.638 -20.732 -21.424 -21.669 -21.764 -21.131 -21.375 -21.172 
249.00 6.2501e-i-43 4.8175e-i-43 -20.783 -19.744 -20.100 -21.052 -21.402 -21.641 -20.986 -21.474 -21.250 
273.90 4.7246e-i-43 3.3740e-i-43 -20.479 -18.700 -19.510 -20.610 -21.069 -21.460 -20.955 -21.445 -21.193 
301.30 3.6861e-i-43 2.4050e-i-43 -20.210 -17.791 -18.959 -20.155 -20.744 -21.227 -21.109 -21.474 -21.290 
331.40 2.4510e-i-43 1.6278e-i-43 -19.766 -17.175 -18.460 -19.670 -20.304 -20.901 -20.665 -20.942 -20.625 
331.40 2.4371e-l-43 1.6453e-l-43 -19.760 -17.516 -18.467 -19.627 -20.293 -20.940 -20.685 -20.864 -20.575 
364.50 1.6297e-i-43 1.1526e-i-43 -19.323 -16.990 -17.944 -19.133 -19.916 -20.659 -20.243 -20.270 -19.774 
401.00 1.1467e-i-43 8.4434e-i-42 -18.942 -16.434 -17.419 -18.704 -19.617 -20.388 -19.953 -19.815 -19.169 
441.10 7.9475e-i-42 5.7846e-i-42 -18.544 -15.721 -16.647 -18.145 -19.241 -20.065 -19.759 -19.545 -18.684 
485.20 5.6956e-i-42 3.8594e-i-42 -18.182 -15.033 -15.792 -17.542 -18.803 -19.720 -19.604 -19.458 -18.379 
645.80 1.0517e-l-42 5.9927e-l-41 -16.348 -12.727 -13.115 -15.085 -16.787 -17.795 -18.048 -17.577 -16.151 
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Table D4. Bolometric and UVOIR luminosities as well as synthetic magnitudes for PISN model HelOOionlNL. 



Age Lboi iuvoiR A^boi Mv Mb My Mr Mi Mj Mh Mk 

(d) (ergs-i) (ergs-i) (mag) (mag) (mag) (mag) (mag) (mag) (mag) (mag) (mag) 

11.55 3.8088e+41 2.4709e+41 -15.245 -14.431 -14.649 -15.228 -15.529 -15.811 -16.213 -16.496 -16.632 
12.71 4.0804e+41 2.6687e-H41 -15.320 -14.324 -14.760 -15.332 -15.631 -15.902 -16.317 -16.611 -16.750 
13.98 4.4089e-H41 2.8803e-H41 -15.404 -14.359 -14.834 -15.418 -15.717 -15.996 -16.415 -16.720 -16.850 
15.38 4.7409e-H41 3.0781e-H41 -15.483 -14.392 -14.887 -15.495 -15.798 -16.083 -16.514 -16.827 -16.942 
16.92 5.0436e-i-41 3.2326e-H41 -15.550 -14.355 -14.917 -15.556 -15.870 -16.157 -16.613 -16.935 -17.036 
18.61 5.3148e+41 3.3305e-H41 -15.607 -14.207 -14.922 -15.598 -15.931 -16.218 -16.716 -17.051 -17.142 
20.47 5.1565e+41 3.0544e-H41 -15.574 -13.646 -14.761 -15.505 -15.886 -16.202 -16.774 -17.130 -17.225 
22.52 4.6197e-H41 2.4974e-H41 -15.455 -12.646 -14.403 -15.264 -15.719 -16.098 -16.766 -17.151 -17.268 

24.77 4.2201e-H41 2.0684e-H41 -15.356 -11.833 -14.031 -15.005 -15.546 -16.000 -16.760 -17.167 -17.315 
27.25 4.1484e-H41 1.8889e-H41 -15.338 -11.393 -13.795 -14.849 -15.463 -15.975 -16.800 -17.217 -17.393 

29.98 4.5362e-H41 2.0219e-H41 -15.435 -11.333 -13.814 -14.903 -15.543 -16.068 -16.913 -17.331 -17.518 

32.99 5.5622e-H41 2.5671e-H41 -15.656 -11.642 -14.116 -15.193 -15.801 -16.291 -17.106 -17.521 -17.701 

36.29 7.4467e-i-41 3.7020e-H41 -15.973 -12.327 -14.628 -15.649 -16.185 -16.613 -17.358 -17.764 -17.923 
39.92 1.0462e-H42 5.6987e-H41 -16.342 -13.543 -15.244 -16.165 -16.620 -16.981 -17.634 -18.026 -18.153 
43.91 1.4940e-H42 8.8947e-H41 -16.729 -14.972 -15.857 -16.664 -17.043 -17.342 -17.907 -18.280 -18.376 

48.30 2.1214e-H42 1.3572e-H42 -17.110 -16.011 -16.404 -17.119 -17.433 -17.681 -18.168 -18.524 -18.592 
53.13 2.9774e-H42 2.0116e-H42 -17.478 -16.779 -16.898 -17.533 -17.798 -18.003 -18.419 -18.765 -18.805 
58.44 4.1165e+42 2.8960e-H42 -17.829 -17.396 -17.350 -17.909 -18.137 -18.306 -18.661 -18.999 -19.015 
64.28 5.5664e-i-42 4.0347e-H42 -18.157 -17.873 -17.758 -18.256 -18.457 -18.592 -18.899 -19.216 -19.220 
70.71 7.4557e-H42 5.5287e-H42 -18.474 -18.337 -18.136 -18.572 -18.754 -18.864 -19.130 -19.430 -19.424 

77.78 9.8045e-i-42 7.3869e-H42 -18.772 -18.748 -18.474 -18.860 -19.029 -19.124 -19.359 -19.639 -19.623 

85.56 1.2666e+43 9.6485e-t-42 -19.050 -19.115 -18.777 -19.124 -19.285 -19.373 -19.584 -19.841 -19.818 
94.12 1.6074e+43 1.2345e-H43 -19.308 -19.438 -19.050 -19.368 -19.526 -19.615 -19.803 -20.032 -20.006 
103.50 1.98966-1-43 1.5376e-i-43 -19.540 -19.705 -19.293 -19.589 -19.748 -19.844 -20.006 -20.205 -20.182 
113.90 2.40516-1-43 1.87056-1-43 -19.746 -19.922 -19.512 -19.790 -19.956 -20.056 -20.191 -20.349 -20.340 
125.30 2.81266-1-43 2.19926-1-43 -19.916 -20.065 -19.700 -19.973 -20.142 -20.244 -20.351 -20.473 -20.474 
137.80 3.14656-1-43 2.46996-1-43 -20.038 -20.101 -19.836 -20.128 -20.304 -20.407 -20.474 -20.561 -20.573 
151.60 3.33676-1-43 2.62156-1-43 -20.101 -20.008 -19.887 -20.238 -20.432 -20.545 -20.550 -20.616 -20.637 
166.80 3.28516-1-43 2.56326-1-43 -20.085 -19.741 -19.790 -20.272 -20.501 -20.645 -20.561 -20.651 -20.670 
183.50 2.95216-1-43 2.25406-1-43 -19.968 -19.291 -19.517 -20.188 -20.463 -20.652 -20.504 -20.670 -20.680 
201.90 2.42496-1-43 1.77836-1-43 -19.755 -18.653 -19.089 -19.970 -20.296 -20.540 -20.382 -20.662 -20.669 
222.10 1.85626-1-43 1.27486-1-43 -19.465 -17.772 -18.535 -19.611 -20.001 -20.335 -20.276 -20.625 -20.637 
244.30 1.36656-1-43 8.68876-1-42 -19.132 -16.763 -17.921 -19.141 -19.631 -20.064 -20.171 -20.515 -20.528 
268.70 9.81056-1-42 5.8724e-i-42 -18.772 -15.845 -17.315 -18.618 -19.229 -19.768 -19.980 -20.275 -20.280 
295.60 7.00626-1-42 4.00586-1-42 -18.407 -14.530 -16.574 -18.024 -18.858 -19.512 -19.756 -19.943 -19.930 
325.20 5.15796-1-42 2.9173e-i-42 -18.074 -13.678 -15.976 -17.526 -18.517 -19.266 -19.536 -19.565 -19.545 
357.70 3.82066-1-42 2.1646e-i-42 -17.748 -13.057 -15.354 -17.030 -18.180 -19.033 -19.311 -19.220 -19.095 
393.50 2.78606-1-42 1.5224e-i-42 -17.406 -12.583 -14.638 -16.490 -17.786 -18.740 -19.044 -18.961 -18.641 
432.90 1.97166-1-42 9.62986-1-41 -17.030 -12.173 -13.989 -15.826 -17.265 -18.306 -18.727 -18.776 -18.224 
476.20 1.34806-1-42 5.88656-1-41 -16.617 -11.785 -13.451 -15.171 -16.724 -17.770 -18.399 -18.288 -17.846 
523.80 8.31696-1-41 3.68546-1-41 -16.093 -11.441 -13.000 -14.608 -16.224 -17.256 -17.869 -17.442 -17.066 
576.20 4.93896-1-41 2.21916-1-41 -15.527 -11.048 -12.511 -14.036 -15.685 -16.693 -17.273 -16.518 -16.173 
633.80 2.77776-1-41 1.22316-1-41 -14.902 -10.618 -11.972 -13.398 -15.044 -16.029 -16.626 -15.610 -15.245 
697.20 1.46686-1-41 6.06786-1-40 -14.209 -10.131 -11.369 -12.681 -14.274 -15.251 -15.931 -14.773 -14.267 
766.90 7.48436-1-40 2.72236-1-40 -13.479 -9.560 -10.688 -11.886 -13.385 -14.359 -15.217 -14.035 -13.202 
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